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13.  ASSTKACT  (M^jumum  200  wcf^) 


In  this  project  we  have  performed  both  experimental  and  theoretical  studies  of 
optical  beam  phase-conjugation  and  of  electromagnetic  scattering  and  propagation  with 
intense  optical  fields.  We  have;  (1)  made  measurements  of  high-power  four-wave  mixing  in 
many  new  polymers  synthesized  at  USC  to  characterize  their  nonlinear  optical  properties; 
(2)  made  what  we  believe  is  the  first  observation  of  nonexponential  attenuation  of  a  weak 
optical  beam  in  a  homogeneous  (highly  scattering)  medium;  (3)  achieved  the  first  time-of- 
flight  measurements  of  electron  drift  velocities  in  photorefractive  insulating  crystals, 
finding  trap-limited  room  temperature  mobilities  -  0.2  cm^  V^s'^  and  conduction  band 
mobilities  over  an  order  of  magnitude  larger,  also  at  room  temperature;  and  (4)  achieved  the 
first  low  energy  (~  picojoule)  optical  logic  gates  in  resonant  atomic  vapors. 
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1.  INTRODUCTION  AND  PROJECT  OBJECTIVES 


Nonlinear  optics,  the  study  of  matter  interactions  with  intense  electromagnetic  fields, 
has  uncovered  surprising  and  useful  physical  effects,  such  as  stimulated  scattering  of  light  by 
which  an  intense  monochromatic  beam  can  be  converted  to  intense  coherent  beams  of 
different  wavelengths.  Recently,  interesting  processes  involving  nonlinear  optical  image- 
bearing  beams  have  been  discovered,  such  as  optical— beam  phase  conjugation.  Optical— beam 
phase  conjugation  is  the  name  given  to  any  process  which  generates,  in  real  time,  the 
time— reversed  replica  of  a  complex,  image— bearing,  optical  beam,  or  which  generates  other 
related  reflected  beams,  both  monochromatic  and  polychromatic.  Phase  conjugation  by 
optical  four— wave  mixing  was  conceived  earlier  in  this  project,  as  were  such  applications  of 
phase-conjugation  as  the  correction  of  image  aberrations,  brightness  enhancement  of  laser 
outputs,  edge-enhancement  of  images,  automatic  steering  of  beams,  and  optical  spectroscopy 
by  Raman-induced  phase  conjugation  (RIPC).  Such  nonlinear  optical  techniques  are 
currently  being  extended  to  perform  many  forms  of  beam  processing,  sorting,  and  routing  on 
picosecond  time  scales.  This  project  aims  at  exploring  and  developing  these  and  other  new 
wave-mixing  processes,  and  the  necessary  nonlinear  materials.  Other  novel  electromagnetic 
scattering  processes  which  occur  at  the  intense  optical  fields  available  from  lasers  are  also 
explored  for  thdr  scientific  and  device  interest.  The  approaches  taken  in  this  project  are  both 
experimental  and  theoretical. 
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2.  MAJOR  ACCOMPLISHMENTS:  12/1/90  to  11/30/91 


The  major  accomplishments  of  this  grant  period  were  as  follows. 


2.1  Synthesis  and  Characterization  of  Optical  Polymers 


About  four  years  ago  in  this  AFOSR  project,  we  began  an  effort  to  develop  rugged  and 
stable  polymer  films,  both  amorphous  and  crystalline,  for  nonlinear  optical  imaging  or 
"real-time”  holography.  This  effort  is  in  collaboration  with  professor  L.  Dalton  and  his 
collaborators  in  the  USC  Chemistry  Deparment  who  synthesize  the  materials  and  perform 
NMR  and  EPR  studies  on  them.  Our  role  in  this  project  is  to  measure:  (1)  one—  and 
two— photon  absorption  spectra,  (2)  film  thickness  using  a  step— profiler  and  electron 
microscopy,  (3)  the  linear  indices  of  refraction  by  reflectivity  measurements,  (4)  space 
components  of  that  part  of  the  degenerate  four— wave  mixing  susceptibility 

tensor  that  respond  in  a  time  less  than  10  ps  (our  instrument  resolution  time)  and  to  measure 
these  vs  u,  (5)  the  parameters  which  determine  the  slower— responding  density  and 
temperature  gratings  and  the  transient  refractive  index  gratings  which  they  create  and  (6)  the 
third— harmonic  light  when  it  is  visible.  To  date  we  have  performed  such  measurements  on 
films  of  several  dozens  of  different  polymers  formed  on  glass  microscope  slides.  We  have 
reported  results  on  over  a  dozen  of  the  more  interesting  polymers  in  various  publications  in 
Sec.  3  and  in  the  conference  presentations  (with  published  summaries)  in  Sec.  4.  In  the 
current  period  we  have  reported  on  refinements  in  the  methods  of  synthesis  and 

□ 

characterization,  details  of  which  are  in  Pubs  3.55  —  3.57.  Pub.  3.55  is  attached  as  Appendix  I . . 
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and  Pub  3.56  is  attached  as  Appendix  II. 
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A  major  accomplishment  of  our  picosecond  polymer  studies  has  been  the  development 
of  straightforward  techniques  for  using  picosecond  pulse  measurements  of  nonlinear  material 
parameters  that  (1)  allow  meaningful  comparison  among  results  from  different  laboratories, 
(2)  allow  realistic  estimates  of  uncertainties  in  the  derived  nonlinear  parameters  that  arise 
from  uncertainties  in  the  transverse  and  longitudional  pulse  profiles,  (3)  allow  full  use  of  the 
sixteen  independent  combinations  of  the  states  of  polarization  of  the  four  beams,  and  (4)  allow 
unambiguous  separation  of  the  various  timescales  of  the  "slow"  responding  nonlinearities  and 
the  relative  contribution  of  the  "fast"  (compared  to  pulse  length)  component.  These 
techniques,  used  in  Appendix  I,  are  described  at  length  in  Pub.  3.44  in  which  we  also  note 
ambiguities  in  the  previous  literature  which  our  technique  can  avoid. 

Other  accomplishments  in  polymer  studies  in  the  current  period  include  succesful 
application  of  the  phase  conjugating  Twyman— Green  interferometer  arrangement  to 
measuring  both  the  real  and  imaginary  parts  of  tensor  elements  of  polymers,  as  we 
described  in  Pubs.  3.56  and  3.57,  the  first  of  which  is  included  here  as  Appendix  II.  We  have 
also  been  able  to  measure  nonlinear  tensor  elements  as  a  function  of  wavelength  from  515  to 
720  nm,  as  shown  in  Appendix  II. 


2.2  Nonexponential  Attenuation  in  EBghly  Scattering  Homogenous  Media;  Artfical 
Kerr  Media 

We  have  made  what  we  believe  is  the  first  observation  of  the  nonexponential 
attenuation  of  a  weak  optical  beam  in  a  homogenous  medium,  that  is,  the  first  exception  to 
the  famous  property  of  weak  light  beams  that  was  first  enunciated  by  Pierre  Bouger  262  years 
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ago:  "In  a  medium  of  uniform  transparency  the  light  remaining  in  a  collimated  beam  is  an 
exponential  function  of  the  length  of  its  path  in  the  medium."  We  expected  such  a  violation 
would  occur  in  oui  latex  suspensions  when  the  attenuation  length  due  to  scattering  was  of  the 
order  of  (or  shorter  than)  the  correlation  distance  between  the  scattering  particles.  Our 
detailed  theory  fit  the  attenuation  measurements  well.  Details  are  given  in  Pub  3.58  which  is 
attached  as  Appendix  III. 

Our  studies  of  nonexponential  beam  attenuation  arose  in  the  course  of  our  longstanding 
project  to  study  homogeneous  suspensions  and  mixtures  of  non— absorbing  particles.  These 
media  exhibit  large  optical  nonlinearity  (Kerr  effect)  by  virtue  of  physical  processes  (such  as 
molecular  reorientation  that  also  cause  a  large  amount  of  light  scattering.  There  is  also  the 
possibility  that,  if  scattering  is  strong  enough,  randomly  distributed  localized  light  modes  will 
form,  similarly  as  "Anderson  localization"  of  electronic  states  forms  in  random  media.  The 
nonexponential  beam  attenuation  is  related  to  light  localization  phenomena. 


2.3  Physics  of  Photorefractive  and  Photoconducting  Insulating  Materials 

In  the  current  project  period  we  achieved  what  we  believe  to  be  the  first  direct 
measurement  of  average  drift  velocity  of  photoexdted  carriers  in  a  photorefractive  crystal  (or 
in  any  insulator)  in  a  known  electric  field.  Details  of  these  measurements  and  comparisons 
other  published  measurements  are  described  in  Pubs.  3.59  to  3.61  which  are  attached  as 
Appendices  IV,  V,  and  VI.  We  outline  these  results  below. 

The  drift  velocities  of  photoexcited  holes  and  electrons  in  known  electric  fields  in 
insulators  have  long  eluded  direct  experimental  observation.  Drift  mobilities  and  carrier 
densities  are  generally  too  low  to  produce  an  observable  Hall  effect.  Transient  photoexcited 
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currents  are  observable  but  they  have  complex  time  envelopes,  and  carrier  densities  can  only 

be  guessed.  Indirectly  inferred  values  of  the  photo-excited  electron  mobility  in  single  n— type 

Bij2SiO20  crystals  have  ranged  from  around  10  to  10  cm  V  s  at  room 

temperature.  We  have  devised  a  holographic  method  to  observe  the  movement  of  a  grating  of 

(»v  10  cm  )  photoexdted  electrons  as  they  drift  in  a  known  applied  field.  The  application 

of  this  to  n— BSO  is  described  in  Appendices  IV,  V,  and  VI.  In  Appendix  IV  we  describe  the 

new  technique  and  report  an  electron  mobility  of  0.24  ±  0.07  cm  V  s  at  room 

temperature.  This  value  is  so  small  as  to  suggest  that  the  electrons  spend  considerable  time  in 

shallow  traps  while  drifting  toward  their  eventual  recombination  with  a  deep  trap.  With  this 

in  mind  we  were  able  to  identify  one  portion  of  our  observed  transient  photocurrent  response 

with  shallow— trap-4imited  drift,  as  we  explain  in  Appendices  V  and  VI.  To  pursue  this  idea 

further  we  are  making  direct  time-of— flight  holographic  measurements  of  mobility  over  a 

temperature  range  (270  to  310  K)  that  show  dear  Arrhenius  forms  for  shallow  trap  depths 

near  300  meV.  A  preliminary  account  of  this  work  is  in  the  published  summary  of  Talk  4.113. 

We  have  also  begun  revising  our  time-of— flight  measurement  technique  for  the  electron 

mobility  so  as  to  see  the  pure  band  conduction  mobility  on  a  short  time  scale  after 

photoexcitation.  In  this  case  the  shallow  traps  play  no  role.  In  Talk  4.11  a  preliminary 
2  —1  —1 

estimate  7  ±  3  cm  V  s  was  given  for  this  mobility  in  a  well  characterized  sample  of 
n-BSO. 

We  also  made  the  first  quantitative  measurements  of  the  spatial  harmonic  content  of 
photorefractive  gratings  and  compared  these  measurements  with  computer  solutions  of  the 
standard  model  of  photoexcitation,  drift,  diffusion  and  direct  recombination  of  a  single  carrier 
spedes.  These  experiments  were  carried  out  in  barium  titanite.  Preliminary  results  were 
reported  in  Talk  4.115. 

In  addition,  we  devised  a  novel  grating  technique  to  measure  the  optical  properties  of 
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the  deep  (ionized  and  occupied)  traps  from  which  the  carriers  are  photo— excited.  We  call  this 
technique  "transport-induced— grating  interferometry."  Its  application  to  the  (unidentified) 
traps  in  cubic  Bij^2'^^O20  outlined  in  Talk  4.116.  Similar  mesurements  have  been 

initiated  to  determine  optical  polarizabilities  of  cubic  n— BSO.  A  preliminary  account  of  first 
results  of  htese  was  given  in  Talk  4.117. 


2.4  Atoms  for  logic 

In  the  current  project  period,  we  began  the  studies  we  proposed  of  a  wholly  different 
class  of  nonlinear  materials  for  image  processing:  thin  cells  of  cesium  or  potassium  vapor  with 
special  "buffering"  by  other  vapors  and  by  wall  collisions.  Preliminary  accounts  of  our  first 
successful  vapor  logic  gates  were  given  in  Talks  4.118  to  4.121.  Some  possible  applications  of 
such  nonlinear  processes  to  the  laser  gravity— wave  detector  were  discussed  in  Talk  4.122. 
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SYNTHESIS  AND  CHARACTERIZATION  OF  THIRD  ORDER  NONLINEAR 

OPTICAL  MATERIALS 
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90089; 
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ABSTRACT 


New  polymers  incorporating  a  variety  of  electroactive  moeities  with 
defined  Ji-electron  conjugation  lengths  have  been  synthesized  and  characterized 
by  degenerate  four  wave  mixing  (DFWM)  techniques.  The  values  for 

these  materials  varied  from  10*^2  jq  10-13  gju  cm.  This  work  has  identified 
several  promising  structures  with  nonlinear  optical  activity  including 
organometallic  and  purely  organic  materials.  The  preparation  of  composite 
materials  has  also  permitted  the  measurement  of  x^^Vot  as  a  function  of  the 
electroactive  unit  concentration. 


INTRODUCTION  i 

Nonlinear  optical  (NLO)  properties  of  polymers  is  a  rapidly  growing  i 

research  field  because  of  the  many  practical  applications  of  polymeric  materials 
in  optoelectronic  devices.^  Polyconjugated  materials  like  polydiacetylene, 
polyacetylene  and  polythiophene  have  been  heavily  investigated.^*^  But  the 
inherent  deficiencies  of  these  materials,  such  as  poor  solubility,  poor  thermal 
and  environmental  stability  and  high  optical  absorption  losses,  prevent  their 
use  in  practical  applications.  Therefore,  new  materials  must  be  found.  Both 
theory  and  experiment  have  demonstrated  that  polymers  with  limited 
conjugated  length  can  manifest  NLO  effects  as  large  as  polymers  with  extended 
polyconjugation.^-6  It  has  been  recognized  that  electron  delocalization  and 
charge  transfer  contribute  to  NLO  effects,  and  that  polymers  with  electron  clouds 
of  enhanced  anharmonicity  are  expected  to  give  rise  to  large  NLO  effects. 

Therefore,  it  is  reasonable  to  test  different  materials  with  different  kinds  of 
constituents  for  optical  nonlinearities.  Such  materials  would  include  those 
containing  transition  metals  and  polymers  with  ionic  or  with  charge  separated 

5 

) 

V 


Abbreviations:  Fourier  Transform  Infrared  (FTIR),  Nuclear  Magnetic 

Resonance  (NMR),  Ultraviolet/ visible  (UV/vis),  Gel  Permeate  Chromatography 
(GPC),  Differential  Scanning  Caloriometer  (DSC),  Dimethylformamide  (DMF), 
Dimethylsulfoxide  (DMSO). 
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species.  With  this  in  mind,  we  have  extended  our  research  to  a  broader  survey 
of  conjugated  electroactive  units.  Here  we  report  part  of  our  effort  in  the 
synthesis  and  characterization  of  new  materials  exhibiting  large  optical 
nonlinearities. 

EXPERIMENTAL 

Polymer  1.  This  polymer  was  synthesized  through  a  precursor  polymer. ^ 
The  mixture  of  prepolymer  (0.904  g,  1.97  mmol)  with  nitrobenzene  (20  mL)  and 
benzoyl  chloride  (3  mL)  was  refluxed  for  4  hrs.  The  reaction  mixture  was 
poured  into  methanol  (50  mL)  and  the  red  precipitate  was  collected  by  filtration 
and  washed  with  methanol  until  the  filtrate  becomes  colorless,  yielding 
polymer  1  (1.09  g,  83%).  NMR  (DMSO-d^)  d  1.6,  (broadened,  6H),  3.9, 
(broadened,  4H),  6.5-8.2  (broadened,  18H)  ppm. 

Anal.  Calcd  for  (C37H28N2Cl206)n:  C,  66.56,  H,  4.20,  N,  4.20,  Cl,  10.64%. 
Found:  C,  66.12,  H,  3.87,  N,  3.75,  Cl,  11.45%. 

Polymer  2.  The  mixture  of  monomer  (l,6-di(l-ethyl,2-methyl- 
quinolidine-6-oxy)  pentane  diiodide)  (l.OOg,  1.36  mmol)  with  triethyl 
orthofomate  (0.739  g,  5.44  mmol)  and  pyridine  (15  mL)  was  refluxed  for  one  hr. 
The  mixture  was  poured  into  THF.  The  precipitate  was  collected  and  washed 
with  cold  methanol.  The  final  product  was  dried  at  40°C  under  vacuum 
overnight,  yielding  polymer  2  (0.73g,  85%).  NMR  (DMSO-d6)  d  1.4,  (t,  J=9.0 
Hz,  CH3-  of  ethyl),  1.8,  (broadened,  -(CH2)3-),  3  0,  (s,  overlapped  with  water  side 
band,  2-methyl),  4.3,  (broadened,  -CXH2-),  5.0  0?roadened,  -CH2-  of  ethyl),  6.5 
(broadened,  =CH-),  7.0-9.0,  (multiple,  aromatic  protons). 

Anal.  Calcd.  for  (Cn9Hi45N808r5)5H20:  C,  56.26,  H,  6.11,  N,  4.41%.  Found: 
C,  56.43,  H,  5.78,  N,  4.40%. 

Polymer  3.  The  mixture  of  monomer^  ^  (1.00  g)  and  AIBN  (0.003g)  was 
added  to  toluene  (20  mL)  in  a  tube.  The  solution  was  frozen  and  thawed  twice 
vmder  vacuum  to  remove  oxygen  and  was  then  sealed.  The  sealed  tube  was 
immersed  in  a  oil  bath  at  70‘’C  overnight.  The  reaction  solution  was  poured 
into  ether  and  the  polymer  was  collected  by  filtration  and  washed  with  ether. 
The  product  was  then  dried  at  room  temperature  under  vacuum  overnight, 
yielding  polymer  3.(0.3  g,  30%). 

Anal.  Calcd.  for  (C3oH32N403)0.3H20:  C,  71.79,  H,  6.50,  N,  11.17%.  Found: 
C,  71.94,  H,  6.54,  N,  11.22%. 

Polymer  4.  Squaric  acid  (0.91  g,  7.98  mmol)  was  dissolved  in  20  mL  hot 
butanol  and  then  N-methylpyrrole  (0.65  g,  7.98  mmol)  was  added.  The  mixture 
was  refluxed  for  2  hrs.  The  dark  blue  precipitate  was  collected  and  washed  with 
water  and  ethanol,  dried  at  50°C  under  vacuum  yielding  polymer  4  (1.16,  80%). 

Anal.  Calcd.  for  (C55H4oN70]o)5H20:  C,  63.03,  H,  4.68,  N,  9.36%.  Found: 
C,  63.07,  H,  4.87,  N,  9.65%. 

Polymer  characterization  and  DFWM  experiments  have  been  described  in 
previous  works. 

RESULTS  AND  DISCUSSION 

The  following  materials  (polymers  and  small  molecules)  have  been 
synthesized  and  characterized  by  DFWM  techniques  (See  Fig.  1). 
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Fig.  1:  Structures  of  NLO  active  materials. 
Polymer  1  is  synthesized  following  reaction  scheme  1. 
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Scheme  1:  Synthesis  of  polymer  1. 


Polymer  1 


The  polymer  structure  has  been  confirmed  by  FTIR,  NMR,  NMR, 
UV/vis  spectroscopic  studies  and  elemental  analysis.  Fig.  2a  shows  the  Uv/vis 
spectrum  of  polymer  1.  Sharp  absorptions  at  ca.  479,  515  and  567nm  can  be 
found.  Polymer  1  is  fully  soluble  in  DMF,  DMSO  and  1-chloronaphthene  and 
optical  quality  films  can  be  cast  from  these  solutions.  GPC  measurements  in 
DMF  show  that  the  weight  average  molecular  weight,  <Mw>/  is  about  20,000 
using  polystyrene  as  a  standard.  This  polymer  is  thermally  stable  up  to  480°C 

Polymer  2  is  synthesized  through  a  well  known  cyanine  dye  formation 
reaction  (see  scheme  2). 


Et 


O. 

^CsHio 


-oa 


CH(OEt), 


NPCH3 

1 1  PyndincorCJumoline 

Et 


Polymer  2 


610 


« 


Both  pyridine  and  quinoline  can  be  used  as  solvents.  Since  these  solvents 
are  not  good  solvents  for  the  final  products,  the  materials  we  obtained  are  not 
high  molecular  weight  polymers.  Elemental  analysis  showed  that  they  have 
only  about  three  repeat  units.  They  are  partially  soluble  in  DMF,  DMSO  and  m- 
cresol.  The  films  can  be  cast  from  solution  or  from  mixed  solution  with 
polycarbonate  on  glass  slides.  NMR  spectra  shows  that  a  broad  peak  at  6.5 
ppm  attributed  to  methine  proton  appears  and  the  2-methyl  peak  of  the 
monomer  decreased  dramatically.  FTIR  spectra  show  similar  features  with 
those  of  model  compound,  6,6'-dimethoxy-quinolidine  cyanine.  Fig.  2b  shows 
the  Uv/vis  spectrum  of  polymer  2  in  which  absorptions  appear  at  540,  580  and 
640  nm. 
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Fig.  2:  Uv/vis  spectra  o^  a).  Polymer  1,  b).  Polymer  2,  c).  Polymer  3,  d). 
Polymer  4  and  e).  Tetraazoanuulene  nickel  complex. 


Polymer  3  is  synthesized  by  a  free  radical  polymerization  of  the 
corresponding  a-vinyl  monomer.!  ^  The  final  polymer  has  MW  of  4000  to  10000 
(polystyrene  as  standard).  NMR,  NMR,  FTIR,  Uv/vis  spectroscopic  data 
and  elemental  emalysis  results  confirm  the  structure  proposed.  Fig.  2c  gives  the 
Uv/vis  sjjectrum  of  polymer  3.  This  polymer  is  thermally  stable  up  to  300®C  as 
shown  by  TGA  studies.  DSC  shows  Tg  at  40°C  and  no  liquid  crystal  transitions 
have  been  observed  between  0®C-200°C. 

Polymer  4  is  a  known  polymer  which  was  synthesized  by  Treibs  and 
Jacobs  with  only  three  repeat  units.  We  modified  the  synthesis  procedure  and 
obtained  a  polymer  with  a  repeat  units  larger  than  five.  The  elemental  analysis 
results  suggest  that  each  squaric  add  unit  is  associated  with  one  water  molecule. 
TGA  and  DSC  studies  support  this  assumption.  A  weight  loss  of  7% 
accompanied  by  a  endothermic  process  were  observed  between  50-140°C.  It  is 
partially  soluble  in  DMF  and  m-cresol  and  forms  blue  solutions.  FTIR  spectra 
show  absorption  at  1775  cm-!(s)  due  to  the  carbonyl  group  of  the  squaric  acid 
unit,  absorption  at  1640  cm*!  (vs)  due  to  the  vinyl  group  in  the  squaric  acid  unit 
and  methyl  absorption  is  observed  around  2890  cm*!.  Fig.  2d  is  the  Uv/vis 
spectrum  of  this  polymer  in  which  the  lowest  absorption  band  appears  around 
610  nm.  This  polymer  is  very  interesting  because  it  can  have  many  resonance 
structures  in  which  the  charge  separated  structure  is  the  most  likely  one.® 
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Tetraazoannulene  complexes  are  synthesized  according  to  a  known 
procedures.^  We  are  interested  in  these  compounds  because  they  contain 
transition  metals  (Ni  can  be  replaced  by  Cu  and  Fe  etc.)  and  they  are  soluble  in 
common  organic  solvents.  They  can  be  derivatized  by  different  substitutents 
both  at  the  R  and  X  positions,  which  provide  opportunity  for  incorporating 
them  into  polymer  backbone.  Fig.  2e  shows  the  Uv/vis  spectrum  of  compound 
with  R=X=H,  in  which  a  window  around  530  nm  can  been  seen.  This  is  an 
advantage  for  DFWM  measurements,  usually  performed  at  532  nm. 

All  these  materials  have  been  demonstrated  NLO  active  by  DFWM 
experiments.  Table  1  summarized  the  results  for  these  materials. 

Table  I.  NLO  properties  of  Polymers  1-4  and  Tetraazoannulene  Nickel 
complex. 

Materials  (esu.  cm)  x  10^3  at  1=532  nm 

polymer  1  0.37 

polymer  2  9.00* 

polymer  2  in  polycarbonate  10%  0.86 

polymer  3  8.50* 

polymer  4  in  polycarbonate  10%  0.90 

complex  in  polycarbonate  10%  2.45 

*  Permanent  grating  has  been  observed  in  high  pulse  energy. 

Fig.  3  shows  a  typical  plot  of  the  DFWM  signal  pulse  energy  versus  the 
delay  time  of  the  backward  pump  beam  B  for  polymer  1.  From  Fig.  3,  we  can  see 
an  extremely  fast  NLO  process  at  zero  delay  time,  which  might  be  attributed  to 
an  electronic  excitation  contribution,  and  a  slower  process  which  might  arise 
from  thermal  contributions.  Polymers  2  and  4  show  similar  signal  patterns.  It 

can  be  noted  from  Table  I.  that  polymer  2  and  polymer  4  have  similar  x^^Va 
values  in  polymer  matrix  which  they  are  larger  than  polymer  1  and,  in  this  case, 
might  mean  that  the  charge  species  and/or  longer  conjugation  make  significant 
contributions.  The  results  in  Table  I  are  comparable  with  those  of 
polyconjugated  polymer  systems  such  as  polyacetylene  and  polythiophene.  It 
demonstrates  that  high  third  order  NLO  effect  can  be  observed  in  polymer 
systems  with  moderate  conjugation  lengths.  In  Table  I,  the  x^^Vot  value  for 
tetraazoanuulene  nickel  complex  is  very  high  and  is  comparable  to  those  of 
phthalocyanine  complexes.^O  We  are  now  trying  to  incorporate  this  complex 
into  a  polymer  backbone  so  that  the  processibility  of  the  materials  can  be 
improved. 

CONCLUSION 

The  processibilities  of  NLO  polymers  can  be  improved  by  introduction  of 
flexible  chain  to  link  electroactive  unit.  The  polymers  so  obtained  with 
moderate  conjugation  lengths  can  manifest  optical  nonlinearities  as  large  as 
polymers  with  extended  polyconjugations. 
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Fig.  3:  Plot  of  DFWM  signal  as  a  function  of  delay  time. 
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ABSTRACT 

Third  order  nonlinear  optical  properties  of  organic  ladder  copolymer  (POL)  system  is  studied 
using  degenerate  four-wave  mixing  with  picosecond  laser  pulse.  Both  the  real  and  imaginary  part 
of  the  third  order  nonlinear  susceptibility  were  determined  by  a  new  phase  conjugate 
interferometric  method  over  the  wavelength  range  of  532  -  720  nm.  From  the  space  symmetry  and 
wavelength  dependence  of  x^^^  we  attribute  the  observed  nonlinearity  to  the  nonlinear 
photoexcitation  of  bipolarion  states  in  this  ladder  copolymer  system. 

1.  INTRODUCTION 


Although  it  is  recognized  that  delocalized  pi-electron  polymers  exhibit  large  nonlinear 
optical  response,  of  particular  interest  is  their  capability  of  exhibiting  fast  nonresonant  x^^^  values 
for  potential  all-optical  device  application.^  However,  several  practical  problems  have  limited  the 
further  exploitation  of  these  nonlinear  optical  properties.  One  of  the  major  difficulty  is  to  find 
large  optical  nonlinearity  with  small  optical  absorption.  True  non-resonant  nonlinearity  is  not 
expect  to  be  large.  Another  avenue  for  exploiting  large  nonresonant  nonlinearity  is  to  look  for 
subgap  resonant  processes  with  small  one  photon  absorption.  Unfortunately  for  subgap  resonant 
process  such  as  two-photon  absorption,  the  dominant  contribution  of  the  optical  nonlinearity 
comes  from  nonlinear  absorption.  2  it  has  been  recently  pointed  out  that  nonlinear  absorption  can 
also  induced  the  undesired  thermal  effects  under  intense  Irradiations  just  as  in  the  case  of  linear 
absorption.  3  in  order  to  ultilize  the  advantage  of  the  fast  reponse  nonllinearity  in  organic 
materials,  it  seems  that  efforts  should  be  concentrate  on  looking  for  large,  real,  subgao  resonant. 
fast  response  optical  nonlinearity  in  organic  materials. 

Ladder  polymers  have  excellent  chemical  and  thermal  stabilities  and  have  recently  been 
shown  that  stable  bipolarion  state  s  can  be  induced  in  model  ladder  structure  related  to  PTL  and  POL 
systems.^  It  is  shown  by  electron  nuclear  double  resonance  (ENDOP^  and  electron  spin  echo  (ESE) 
spectroscopies,  that  electron  delocalization  might  not  extended  much  beyond  one  ladder  repeat  unit. 
In  the  view  of  the  fact  that  the  polarionic  electron  delocalization  in  ladder  polymers  may  in  fact  be 
confined  to  an  average  of  one  repeat  unit,^  it  seemed  appropriate  to  study  the  possibility  of 
formation  of  the  less-confined  blpolarionic  states  in  ladder  model  systems  such  as  POL  structure 
Illustrated  in  Figure  1 . 

Experimental  results  of  electro-chemical  doping  have  indicated  the  bipolarion  state  is  the 
more  stable  state  and  can  easily  be  achieved  with  excess  dopant.  In  addition,  electron-donating 
groups  stabilized  the  blpolarionic  state  significantly.  Bipolarion  states  generated  by  electro¬ 
chemical  doping  in  monomer  PTL  ladder  model  system  has  been  shown  to  be  stable  in  the  presence 
of  moist  air  for  several  days.'^ 
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In  this  paper  .  we  study  the  possibility  of  photoexcitation  of  bipolarions  in  POL  model 
ladder  system  by  investigating  optical  nonlinearity  in  this  system.  Both  the  real  and  imaginary 
part  of  the  third  order  nonlinearity  is  determined  by  using  a  modified  phase  conjugate 
interferometric  method  over  the  wavelength  range  of  532  -  720  nm.  comparison  between  the 
nonlinear  photoexcitation  of  bipolarons  and  the  bipolarions  generated  by  electro-chemical  doping 
will  be  discussed. 


2.  EXPERIMENT  AND  RESULTS 

Among  different  experimental  methods  in  studies  of  third  order  nonlinear  optical  properties 
of  materials,  degenerate  four-wave  mixing  is  the  most  sensitive  one.  Extreme  small  index  change  ( 
An  =  10'®  )  can  be  detected  easily  by  this  method.  However  this  method  is  only  sensitive  to  the 
magnitude  of  the  complex  third  order  nonlinear  susceptibilities  >  the  phase  of  cannot  be 
detected  by  this  method.  In  this  section  we  described  a  new  interferometric  method  which  can 
determined  the  real  and  imaginaiy  part  of  (-“•  w.  (o,  -co)  accurately  at  the  same  time.  This 
technique  was  first  proposed  by  Rentzepis  et.al.  ®  In  this  study  we  have  extended  this  technique  to 
the  picosecond  pulse  measurement  in  nonlinear  thin  films. 

This  interferometric  method  is  base  on  a  modified  Twyman-Green  interferometer,  which 
includes  two  phase  conjugate  mirrors.  As  will  be  discussed  in  the  following  ,  this  interferometer 
makes  possible  the  accurate  measurement  of  relative  phase  shift  between  two  phase  conjugators.  A 
schematic  representation  of  this  double  phase  conjugate  interferometer  is  shown  in  Figure  2,  two 
phase  conjugate  mirrors  are  based  on  degenerate  four-wave  mixing  scheme.  The  forward  pump 
beam  and  probe  beam  are  split  into  two  sets  by  a  beam  splitter  BS,  and  follow  the  two  paths  shown  in 
Figure  2.  The  backweird  pump  beams  are  provided  by  the  back  reflections  of  the  forward  pump  beams 
from  the  mirrors  right  next  to  the  two  nonlinear  samples.  After  interacting  at  the  phase  conjugate 
mirror,  phase  conjugate  beams  PCI  and  PC2,  and  two  retroreflected  backward  backward  pump 
beams  recombine  at  BS.  Such  a  setup  consists  two  Twyman-Green  interferometers,  two  backward 
pump  beams  interfere  with  each  other  forms  a  conventional  interferometer  which  we  designate  it  as 
pump  interferometer  PI.  The  two  phase  conjugate  beams  constitute  a  phase-conjugate 
interferometer  Cl.  The  edignment  of  one  interferometer,  PI  or  Cl,  results  in  automatic  alignment  of 
the  other.  In  the  experiment  .  the  alignment  of  PI  is  achieved  by  eliminating  the  the  interference 
fringes  in  front  of  the  detector  PDl,  perfect  alignment  is  reach  when  a  circular  ring  pattern  is  found. 
The  center  spot  of  this  circular  pattern  is  selected  by  a  pinhole  and  then  sent  to  the  detector.  A 
piezoelectric  transducer,  PZT  is  mount  on  the  back  of  the  high  reflector  in  one  arm  to  vary  the  path 
length  difference  between  two  arms.  Also  path  difference  between  two  arms  must  be  roughly  the 
same  in  order  to  temporal  overlap  two  picosecond  interference  pulses.  (In  our  case  pulse  length  is 
about  6mm) 

The  best  way  to  understand  the  data  from  this  modified  phase  conjugate  interferometers  is  to 
introduce  a  variable  phase  shift  A(p  into  one  of  the  arms  of  the  interferometer  and  analyze  how  the 
interference  pattern  change  with  this  phase  shift.  This  variable  phase  shift  can  be  introduced  by 
moving  one  of  the  backreflection  mirrors  behind  the  nonlinear  sample.  The  interference  signal  of 
the  pump  interferometer,  Ipj  is 

Ipi  =Il+  l2  -  2(Ii  I  2)^'^^cos(47t(li  -I2  )A+A(p),  (1) 

where  Ij  and  1  2  are  the  intensities  of  the  returning  pump  beams  that  follow  the  optical  path 
lengths  21 1  and  21  2  of  arms  1  and  2.  respectively.  Using  the  expression  for  phase  conjugate  signal, 
it  is  easy  to  show  that  the  interference  signal  Ici  of  the  phase  conjugate  interferometer  Cl  can  be 
written  as 


^C2'  2(Ici  I  c2)^/2cos(47t(li  -1  2  )/^  +  -  <P02l  • 


(2) 
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where  (poi  and  (po2  are  absolute  phase  shifts  of  the  phase  conjugators.  which  depend  on  the  type  of 
diffraction  grating  created  by  the  pump  and  probe  beams.  Comparison  of  equations  (1)  and  (2)  shows 
that  the  interference  signals  at  two  outputs  are  described  by  the  same  function  of  the  induced  phase 
shift  A<j).  The  only  difference  is  the  phase  shift  <pqi  -  <po2-  'This  phase  shift  is  equal  to  the  difference 

of  absolute  phase  of  (-o).  co.  (o,  -w)  in  the  two  nonlinear  mediums.  This  means  that,  if  the  absolute 
phase  shift  in  one  of  the  nonlinear  material  is  known,  then  the  absolute  phase  of  other  materials 
can  be  determined  easily.  In  the  case  when  same  medium  is  used  in  both  arms,  the  interference 
signal  at  both  outputs  should  be  exactly  in  phase.  Also  the  ration  between  the  maximum  and 
minimum  of  the  interference  pattern  can  be  used  to  determine  the  ratio  of  I'W.  (o,  co,  -co)  in  both 
mediums  accurately. 

To  determine  the  absolute  phase  of  in  POL  ladder  copolymer  system,  experiments  using  this 
interferometer  have  been  performed  over  the  wavelength  range  of  532  -  720  nm.  To  check  the 
interferometer  setup  we  initially  used  CS2  in  both  interferometer  arms.  Two  1mm  thick  cells  filled 
with  CS2  were  placed  in  two  arms.  In  the  Figure  3,  the  interference  signal  of  both  Cl  and  PI  are 
plotted  as  a  function  of  the  PZT  voltage.  These  data  shows  two  signals  are  exactly  in  phase.  ( cpoi  - 
cp()2  =  ±  3  ®  )  Next  we  put  POL  ladder  copolymer  thin  fUms  in  one  arm,  and  CS2  cell  in  the  other.  CS2 
is  transparent  over  this  wavelength  range,  and  it  is  well  understood  that  its  nonlinearity  is  mostly 
real  and  positive.  Experimental  results  at  590  nm  are  shown  in  Figure  4,  from  this  figure  we  found 

two  interference  signals  are  shifted  by  114  ±  3*^  .  this  suggests  that  x'^^  ("W-  "“)  of  POL 

copolymer  at  590  nm  is  mostly  imaginary,  and  have  a  negative  real  part. 

To  determine  the  magnitude  of  the  complex  x^^^  'W)  .  we  have  performed  the 

traditional  four-wave  mixing  experiment  in  POL  ladder  copolymer  system.  In  this  case,  the 
backward  pump  beams  is  provided  by  a  third  laser  pulse  instead  of  using  the  reflection  from  the 
mirrors.  Knowing  the  magnitude  and  the  phase  of  x^^^  (-o).  to,  to,  -co),  we  can  determine  its  real  and 
imaginary  part  at  the  same  time.  Also,  by  this  setup,  we  are  able  to  time  delay  pulse  of  different 
beams  to  study  the  transient  dynamics  of  the  photoexcited  bipolarionlc  states.  A  laical  decay  data 
for  phase  conjugate  signal  as  a  function  backward  pump  beam  delay  is  shown  in  Figure  5. 

We  also  study  the  tensor  symmetry  of  the  X^^^jkl  studying  the  polarization 

dependence  of  the  phase  conjugate  signal.  In  the  case  when  the  probe  beam  polarization  is  vertical 
while  the  pump  beam  polarizations  are  horizontal,  phase  conjugation  signal  is  at  least  1000  times 
smaller  than  in  the  case  when  all  beam  polarization  are  horizontal.  This  result  indicates  that 
X^^^1221  X^^^ll  11/30  ,  therefore  suggest  that  the  observed  nonlinearity  is  relate  to  a  nonlinear 
population  excitation. 


3.  DISCUSSION 


In  order  to  understand  the  the  photogeneration  bipolaron  process,  we  will  first  discuss  some 
feature  of  bipolaron  generation  by  electro-chemical  doping.  In  essence  ,  during  the  doping  prcx;ess, 
a  new  subgap  band  structure  which  is  shifted  far  to  the  red  ,  replace  the  original  7t-7t*  band  of  the  POL 
compounds.  There  are  three  significant  signatures  of  the  generation  of  bipolarionlc  states:  (1)  The 
formation  of  localized  structural  distortions  with  associated  localized  vibrational  modes  (IRAV 
modes)  in  the  midinfrared.  (2)  The  generation  of  symmetric  gap  states  and  the  associated  electronic 
transitions.  (3)  The  reversed  spin-charge  relation,  i.e.,  charge  storage  in  spinless  blpolarons. 

In  POL  ladder  system  all  these  features  have  been  verified  experimentally  during  electro¬ 
chemical  doping.  Absorption  spectra  change  due  to  doping  indicates  the  formation  of  subgap  states. 
^  We  have  also  performed  electron  spin  resonance  spectroscopy  (ESFd  experiment  for  the  doped  and 
undoped  samples,  no  ESR  signal  was  observed  under  the  saturate  doping  case.  All  these  results  in 
doped  POL  systems  are  in  consistent  with  the  picture  of  blpolarons  generation. 
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Next  we  discuss  the  process  of  photoexcitation  of  bipolaronic  state.  The  direct 
photogeneration  of  charge  bipolaron  is  clearly  Impossible.  After  electron-hole  photoinjection,  the 
confinement  resulting  from  the  nondegenrate  ground  state  wOl  lead  to  rapid  formation  of  a  neutral 
"bipolaron",  l.e,,  a  neutral  exclton  with  one  electron  in  the  lower  gap  state  and  one  electron  in  the 
upper  state.  Such  neutral  "blpolarons"  can  be  expected  to  undergo  rapid  radiative  decay.  This 
probably  is  the  dominant  mechanism  at  earlier  times.  Because  of  the  transverse  bandwidth  due  to 
interchain  coupling,  some  of  the  initial  photoexcitations  occurs  as  electron  and  holes  on  different 
chains.  Chain  distortions  will  quickly  form  around  these  single  charges,  leading  to  polaron 
formation.  Bipolcirons  can  then  be  generated  if  energetically  favorable. 

From  the  results  of  our  studies  of  optical  nonlinearity,  first  of  all,  the  tensor  symmetry  { 
X^^^1221  <  111/30  )  of  the  observed  nonlinearity  indicated  that  the  nonlinear  mechanism  is 

relate  to  photoinduced  population  of  certain  "excitons".  This  symmetry  holds  for  all  the  wavelength 
we  studied. 

In  Figure  6  ,  we  plot  the  real  and  imaginary  part  of  the  /  a  as  a  function  of 

wavelength.  The  solid  curve  is  deduced  from  the  absorption  spectra  change  during  the  doping.  By 
expressing  the  results  this  way,  we  have  normalized  the  observed  nonlinearity  to  the  photoinduced 
change  of  optical  properties  (refractive  index  or  absorption)  per  photon  absorbed  by  the  POL  ladder 
copolymer  thin  films.  This  quantities  is  independent  of  the  POL  molecular  density  in  the  sample 
and  can  be  compare  with  other  sample  with  different  POL  concentration  or  even  monomer  unit 
mixed  in  polymer  matrix.  As  seen  in  this  figure  ,  the  Imaginary  part  of  the  nonlinearity  ,  or  in  other 
word  ,  the  photoinduced  absorption,  agrees  well  with  the  wavelength  dependence  of  the  absorption 
change  due  to  electro-chemical  dof  ing.  This  indicates  the  photoinduced  absorption  spectrum  is 
similar  to  the  doping-induced  absorption  spectrum,  therefore  indicates  the  photogeneration  of 
bipolaronic  state  in  POL  ladder  copolymers. 

We  have  also  compared  the  results  of  copolymer  and  dilute  monomer  solid  solution,  in  the 
case  of  monomer  solution,  about  3  times  smaller  x'^^  (-O),  o),  (o,  -to)  /  a  value  is  found.  This  indicates 
the  contribution  of  interchain  coupling  in  the  high  density  copolymer  structure.  However,  5  times 
higher  optical  damage  threshold  (2GW/cm2  for  5  psec  pulse)  was  found  in  the  monomer  solid 
solution  (0.5%  mol/mol).  Even  in  this  dilute  concentration  of  POL  units  ,  no  evidence  of  polaron 
generation  is  observed. 

Transient  decay  dynamics  of  the  photogeneration  of  bipolarons  is  studied  by  the  picosecond 
resolved  (5  psec  resolution)  degenerate  four-wave  mixing  experiment.  As  seen  form  Figure  5,  two 
decay  processes  is  evident.  A  pump  intensity  dependent  fast  decay  (  x  =  8  ±  2  psec  at  intensity  of  0.2 
GW/cm^),  and  an  intensity  independent  slow  decay  (  t  =120  ±  5  psec).  The  intensity  dependent  fast 
decay  characterized  the  rapid  recombination  of  the  neutral  bipolarons.  The  slow  decay  process  may 
be  relate  to  interactior.  of  defects  in  the  structure.  Same  type  of  delay  behavior  was  observed  when 
the  forward  pump  pulse  was  delayed.  This  result  indicate  s  that  the  decay  processes  are  local 
interaction  coupling. 

The  significance  of  this  observation  is  underlined  by  the  fact  that  effective  photogeneration 
of  bipolaronic  states  can  be  achieved  in  short  conjugated  ladder  systems  such  as  POL.  Our  results 
Indicate  the  possibility  of  obtaining  large  subgap  resonant  "bipolaron  enhanced  x^^^  value  in  one 
repeat  unit  of  ladder  polymer  structure.  And  as  seen  in  Figure  6  ,  unlike  two  photon  absorption  , 
large  real  contribution  of  subgap  resonant  nonlinearity  can  be  obtain  at  certain  wavelength 
windows. 

Because  of  the  fact  that  effective  nonlinear  excitation  process  can  be  support  in  short  ladder 
segments,  we  are  able  to  observed  high  nonlinear  optical  activity  in  copol)mier  structures.  The 
advantage  of  using  copolymer  structure  is  the  nonactive  spacer  Improved  the  solubility  of  the 
pol)nner.  Thus  processiblllty  is  attained  while  achieving  a  controllable  optical  widow  and  high 
nonlinear  optical  activity.  We  are  currently  exploring  options  of  controlling  optical  windows,  film 
forming  chai  acteristic  and  thermal  properties  of  different  ladder  copolymer  systems. 
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4.  CONCLUSION 


In  summary,  we  have  observed  the  first  experiment  evidence  of  photogeneration  of 
bipolaronic  states  in  organic  ladder  copolymer  systems.  The  photoinduced  absorption  spectrums  is 
similar  to  the  doping-  induced  absorption  spectrums.  Large,  supgap  resonant,  real,  fast  response 
optical  nonlinearity  together  with  good  thin  film  processiblllty  and  stability  can  be  achieved  in 
these  copolymer  systems. 
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'  lirror 


Interference  Signal  (a.u.) 


Fig.  3  Interference  patterns  of  the  pump  interferometer  (—#--) 
and  phase  conjugate  interferometer  ( ---A—  ) 
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Interference  Signal  (a.u.) 


Piezo  Voltage  (V) 

Fig.  4  Interference  patterns  of  the  pump  interferometer 
(--#—)  and  the  phase  conjugate  interferometer  (---A--). 
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Phase  conjugate  signal  (a.u.) 


Delay  (psec) 

Fig.  5  Phase  conjugate  signal  as  a  function  of  the  delay  of 
backward  pump  pulse. 
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Fig.  6  Wavelength  dependence  of  the  real  (—A—)  and  Imaginary  (-— o-— )  part 

of  jr(3)  (-0),  (0,  CO,  -CO)  /  a  in  POL  ladder  copolymer  1.  The  solid  line  is  the 
normalized  absorption  spectra  change  during  electro-chemical  doping. 


f24  /  SPIE  Vol.  1 337  Nonlinear  Optica!  Properties  of  Organic  Materials  m  fl 990) 


PHYSICAL  REVIEW  B 


VOLUME  43,  NUMBER  16 


TL 

1  JUNE  1991 


Anomalous  optical  transmission  in  homogeneous  latex  suspensions 

Patrick  W.  Tam  and  R.  W.  Hellwarth 

Departments  of  Physics  and  Electrical  Engineering,  University  ofScr’thern  California, 

Los  Angeles,  California  90089-0484 
(Received  15  October  1990) 

We  have  observed  a  deviation  from  the  law  of  exponential  attenuation  of  a  633-nm  optical  beam 
in  a  homogeneous  scattering  medium  in  the  low-intensity  limit,  where  nonlinear  optical  effects  or 
thermal  effects  are  insignificant  (Bouguer's  law).  This  deviation  occurs  when  the  spatial  correlation 
length  is  comparable  to  or  larger  than  the  asymptotic  attenuation  length,  whence  the  assumption  of 
independent  scattering  events  is  violated.  We  develop  a  statistical  model  for  the  optical  dielectric 
fluctuations  which  fits  well  our  nonexponential  beam-transmission  data  in  homogeneous  suspensions 
of  l-/im  latex  spheres.  We  use  this  to  determine  the  correlation  lengths  (7  to  12  pm)  of  the  dielectric 
fluctuations  in  these  suspensions.  These  correlation  lengths  are  consistent  with  values  that  we 
deduce  from  the  observed  angular  spreads  of  forward  light  scattering. 


I.  INTRODUCTION 

In  1729  Pierre  Bouguer  proposed  that  “in  a  medium  of 
uniform  transparency  the  light  remaining  in  a  collimated 
beam  is  an  exponential  function  of  the  length  of  its  path 
in  the  medium,”'  an  idea  which  was  popularized  by  J.  H. 
Lambert  thirty  years  later.  Here  we  report  an  observed 
exception  to  this  law  of  exponential  attenuation  for  low- 
intensity  light  beams.  We  study  liquid  suspensions  in 
which  the  attenuation  length  (  ^  1-3  ^m)  is  of  the  same 
order  of  magnitude  as  the  correlation  length  (^7-12  pm) 
of  the  fluctuations  of  the  optical  dielectric  constant,  thus 
violating  the  usual  assumption  (stated  explicitly  by 
Bouguer)  that  events  which  attenuate  the  beam  are  un¬ 
correlated.  We  measure  the  relation  between  the  optical 
transmission  T  through  latex  suspensions  (polystyrene 
spheres  in  water)  and  their  thickness  L  for  five  concentra¬ 
tions  of  the  same  particle  size  for  which  the  transmission 
does  not  scale  exponentially  with  the  thickness.  Optical 
transmission  in  latex  suspensions  has  been  studied  by 
Ishimaru  and  Kuga,  who  studied  the  relation  between  the 
transmitted  beam  fraction  T  of  samples  of  latex  suspen¬ 
sions  and  the  latex  particle  density.’  They  observed  ex¬ 
ponential  decays  of  transmission  T  with  thickness  L,  but 
the  attenuation  coefficient  was  not  proportional  to  parti¬ 
cle  density  N  at  the  higher  densities.  That  is,  “Beer’s 
law”  was  violated,  a  fact  with  which  our  results  will 
correlate.  The  error  in  attenuation  measurements  arising 
from  scattered  light  in  latex  suspensions  at  low  concen¬ 
trations  (<0.27%)  was  studied  by  Zaccanti  and  Brus- 
caglioni.'  This  error  our  method  will  avoid.  In  this  pa¬ 
per,  we  describe  an  experimental  technique  that  is  sensi¬ 
tive  to  nonexponential  transmission,  and  then  develop  a 
Gaussian  statistical  model  for  dielectric  fluctuations 
which  fits  well  the  nonexponential  attenuation  we  ob¬ 
serve.  We  determine  the  spatial  correlation-length  pa¬ 
rameter  L,  in  this  model  which  gives  the  best  fit  to  the 
data  using  two  different  forms  for  a  correlation  function. 
These  forms  give  L,  values  which  agree  within  20%  with 


each  other  and  with  the  correlation  length  estimated 
from  the  angular  spectrum  of  the  forward  singly- 
scattered  light.  We  describe  numerous  checks  of  possible 
extraneous  influences  and  find  these  negligible.  Our 
correlation-length  measurements  suggest  the  existence  of 
long-range  forces. 

II.  MEASUREMENT  OF  BEAM  TRANSMISSION 

We  use  as  scattering  material  colloidal  suspensions  of 
monodisperse  latex  (i.e.,  polystyrene)  spheres  of  diameter 
1.09±0.03  pm  as  specified  by  Sigma  Chemical  Company. 
As  purchased,  the  volume  fraction  is  approximately  10% 
which  corresponds  to  a  number  density  A  =  1.48X10" 
cm~  \  and  a  nearest -neighbor  separation  of  1.9  pm  (as¬ 
suming  a  simple  cubic  lattice),  which  is  about  twice  the 
diameter  of  the  particles.  Using  a  pH  meter’^  we  mea¬ 
sured  this  sample  to  have  a  pH  of  6.6±0.3.  We  measured 
the  dc  conductivity  to  be  (7.3±0.5)X  10~'*n~'  m“'. 
More  dilute  samples  are  prepared  by  adding  deionized 
water.  More  concentrated  samples  are  prepared  by  let¬ 
ting  the  suspension  settle  and  decanting  the  water.  The 
samples  are  homogenized  in  an  ultrasonic  bath  for 
several  minutes,  after  which  they  remain  stable  for  more 
than  a  half  day,  adequate  time  for  a  complete  set  of  data 
to  be  obtained. 

We  have  devised  an  experimental  technique  for  deter¬ 
mining  beam  transmission  T  vs  sample  thickness  L  that 
employs  a  direct  comparison  with  the  transmission 
through  a  variable-length  cell  containing  a  dye  of  known 
transmission  function.  This  comparison  is  essentially  a 
null  comparison  that  keeps  all  detector  signals  within  a 
small  range.  The  experimental  setup  is  described  in  Fig. 
1.  The  variable-thickness  spectrometer  cell  Cl  contains 
the  latex  suspension  while  the  standard  comparison  cell 
C2  contains  a  calibrated  dye  solution  (Nile  Blue  in 
methanol).'  A  He-Ne  laser  emits  a  vertically  polarized 
10  mW,  633  nm  TEM,^,  Gaussian  beam  of  radius  0.7  mm 
at  I  /c'  intensity  whose  attenuation  in  the  cells  is  studied 
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FIG.  1.  Experimental  setup  for  the  measurement  of  optical 
transmission.  Beam  splitter  BS  and  photodiode  PDl  monitor 
the  power  incident  on  the  variable-length  cell  Cl  containing  the 
latex  suspension.  The  compensating  variable  attenuator  C2 
(containing  dye  solution)  keeps  the  transmitted  beam  power  into 
the  (500  channel)  optical  multichannel  analyzer  OMA  nearly 
constant.  Filter  F  blocks  fluorescence  from  the  dye.  Data  is  an¬ 
alyzed  by  digital  oscilloscope  DS. 


as  follows.  A  small  fraction  of  this  beam  is  reflected  off  a 
beam  splitter  BS  to  a  monitoring  photodiode  PDl,  the 
output  P,  which  is  displayed  on  the  digital  oscillo¬ 
scope  DS.  The  other  part  of  the  beam  traverses  cells  Cl 
and  C2,  after  which  the  beam  (of  power  ~  10  nW)  is 
detected  by  an  optical  multichannel  analyzer  OMA 
placed  two  meters  beyond  the  latex  sample.*  The  OMA 
measures  the  intensity  profile  of  the  transmitted  beam  at 
angles  0  to  the  beam  axis  up  to  several  beam  widths.  The 
OMA  is  shielded  from  fluorescence  by  interference  filter 
F.  The  signal  T)  from  the  OMA  is  averaged  for  25 
seconds  and  then  stored  in  the  digital  oscilloscope  DS  as 
the  intensity  I  opposite  its  address  A',  (which  is  a  measure 
of  the  angular  deviation  6  from  beam  axis).  A  typical 
stored  trace  is  shown  in  Fig.  2  which  shows  the  transmit¬ 
ted  Gaussian  beam  superposed  on  the  background  B  of 
scattered  light.  Here,  as  in  all  of  our  transmitted  beam 
profile  measurements,  we  could  never  see  any  deviation  of 
the  beam  shape  or  width  from  that  which  we  observed 
with  the  cell  Cl  absent. 

By  rotating  the  threaded  outer  shell  of  the  cell  Cl,  we 
translate  one  of  its  glass  windows  through  a  distance 
given  by  the  difference  between  two  micrometer  readings 
(correct  to  within  ±0.25  fim)  on  the  outer  shell  which  are 
indicated  by  S  in  Fig.  1.  We  express  the  relation  between 
the  scale  reading  5  on  Cl  and  the  actual  thickness  L  of 
the  suspension  by  S  =  L  ±L,|,  where  Lj,  is  the  scale  read¬ 
ing  corresponding  to  zero  thickness.  To  avoid  backlash 
problems,  we  take  data  by  decreasing  the  thickness  L  of 
Cl,  while  increasing  the  thickness  Lj  of  C2  (with  corre¬ 
sponding  scale  reading  ±L^(,)  to  keep  the  peak 

signal  at  the  OMA  center  channel  constant  to  within 
30%.  This  reduces  the  error  in  the  measurements  due  to 
the  nonlinear  response  of  the  detector  to  less  than  1%. 
We  start  from  the  thickness  of  the  suspension  in  Cl 
where  we  can  detect  a  transmitted  signal  and  reduce  L  to 
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FIG.  2.  A  typical  trace  of  beam  intensity  1  vs  angle  0  from 
beam  axis  as  measured  by  the  OMA.  The  background  B  of  scat¬ 
tered  intensity  is  approximately  constant  in  this  small  angular 
range  and  is  computed  by  averaging  the  intensity  from  point  1 
to  2.  The  area  A  under  the  Gaussian  curve,  with  the  back¬ 
ground  B  as  the  baseline,  is  computed  by  a  processor  in  the  os¬ 
cilloscope  and  used  to  calculate  the  beam  fraction  T  transmitted 
through  the  suspension,  as  described  in  the  text  and  plotted  in 
Fig.  3.  This  trace  is  that  used  to  obtain  the  data  point  marked  P 
on  trace  B  in  Fig.  3.  All  other  data  points  showed  the  same 
beam  shape  and  width. 


the  thickness  where  the  transmission  levels  off  with  fur¬ 
ther  rotation  of  the  outer  shell.  After  each  adjustment  of 
the  thickness  L,  the  suspension  is  allowed  to  stabilize  for 
several  minutes  before  the  measurement  is  performed. 
This  is  required  to  achieve  reproducible  data.  At  each  of 
the  scale  settings  5  and  Sj,  a  processor  in  the  oscilloscope 
computes  the  average  background  intensity  B  using  the 
formula  5  =2^ (^)/(number  of  points  sampled).  Here 
the  summation  is  taken  from  the  point  1  to  2  as  shown  in 
Fig.  2.  The  area  A  (5)  under  the  Gaussian  beam  profile  is 
then  computed  by  /4  (5)  =  ]£'[/(0)  — fi]  where  the  sum¬ 
mation  2’  is  taken  from  point  2  to  3. 

We  calibrate  the  dye  solution  by  putting  it  in  cell  C2 
with  Cl  removed.  A  photodiode  is  used  instead  of  the 
OMA  and  the  transmitted  power  is  measured  as  the 
thickness  of  C2  is  decreased,  with  calibrated  neutral  den¬ 
sity  filters  placed  in  front  of  the  detector  when  necessary. 
The  transmitted  power  obeys  Bouguer’s  law  with  an  ab¬ 
sorption  constant  0^=0. 501  ±0.001  mm  '  calculated  by 
a  least-squares  fitting.  To  determine  the  transmission  of 
the  suspension  TiS)  from  and  zf  (S)  it  is  necessary  to 
measure  the  absolute  value  of  this  transmission  at  some 
reference  scale  reading.  Therefore,  at  the  end  of  each  set 
of  measurements  for  a  particular  particle  density,  we 
determine  the  transmission  T{S,„)  at  the  minimum  scale 
reading  S,„  attainable  (with  corresponding  scale  reading 
on  cell  C2)  by  removing  the  cell  C2  and  measuring 
the  transmitted  beam  power  P(S,„  )  with  a  silicon  detec¬ 
tor  (here  the  transmitted  power  is  too  high  for  the  OMA). 
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The  cell  Cl  is  then  removed  and  the  incident  beam  power 
Pq  measured  with  the  same  detector.  Taking  into  ac¬ 
count  4.5%  reflection  loss  R  from  the  window  (4.2% 
from  the  glass-air  interface  and  0.3%  from  the  glass-fluid 
interface),  we  calculate  the  transmission  using 
T{S„)  =  P{S„)/[Pq{\—R)~].  From  this  we  get  the 
transmission  at  other  scale  readings,  to  within  4%,  by  the 
relation  T{S)=T{S„)A  (S)exp[aj(Sj„  —  )]/.4  {S,„ ). 

At  the  end  of  the  experiment  on  a  given  suspension  the 
sample  is  weighed.  After  being  dried  in  a  vacuum  oven 
for  a  day,  the  dried  polystyrene  is  weighed  again  to  deter¬ 
mine  the  weight  fraction  of  the  polystyrene  spheres.  Us¬ 
ing  the  density  of  a  single  latex  sphere  of  1.05  gcm“^ 
quoted  by  Sigma  Chemical  Company,  we  calculate  the 
particle  number  density  N  of  each  sample  to  within  an  es¬ 
timated  uncertainty  of  less  than  13%.  The  results  of  this 
procedure  for  five  densities  large  enough  to  show  devia¬ 
tions  from  Bouguer’s  law  are  shown  in  Fig.  3.  The  devia¬ 
tion  from  exponential  behavior  is  most  evident  for  the 
higher  densities.  The  lines  in  Fig.  3  are  the  best  fit  to  an 
equation  which  we  derive  using  a  phenomenological 
model  for  the  statistics  of  the  spatial  distribution  of  the 
spheres.  This  equation  fits  the  data  to  within  experimen¬ 
tal  error  for  parameters  which  we  describe  in  the  next 
section.  Since  the  cell  has  to  be  disassembled  for  cleaning 
after  each  set  of  measurements,  the  scale  reading  changes 
from  one  set  of  measurements  to  another  by  a  small 
amount.  The  scale  reading  Lq  at  zero  thickness,  which 
changes  from  one  set  of  data  to  the  other,  is  treated  as  a 
fitting  parameter  in  the  model  which  we  now  develop. 


III.  THEORETICAL  MODEL  FOR  BEAM 
TRANSMISSION 


We  assume  that  a  Gaussian  beam  propagates  along  the 
z  axis  and  enters  a  turbid  medium  (of  thickness  L)  which 
can  be  characterized  by  a  real  scalar  dielectric  constant 
e(x,/).  We  also  assume  that  the  optical  electric  vector 
can  be  written  as  Re[x£(x,r  )exp(/kz  — /&>/ )]  where  E 
obeys  the  modified  paraxial  wave  equation: 


— ^  +  ^7+2'^^  -f  —  8e(x,t ) 
“■  3y  ^  3z  c  - 


£(x,/  )  =  0 


(1) 


Here  Se(x,/ )  =  e(x,f )  — is  the  zero-mean  fluctuating 
part  of  the  dielectric  constant  that  arises  because  of  the 
spheres  in  the  water.  The  thickness  L  in  our  experiment 
is  in  the  range  of  10-50  pm  so  that  the  light  sees  a  purely 
static  medium  (no  time  dependence)  on  a  single  pass. 
Also  we  can  neglect  the  transverse  diffraction  of  the  beam 
inside  the  medium.  If,  in  addition,  we  can  neglect  the 
contribution  from  scattered  light  to  the  electric  field  of 
the  beam  at  the  exit  plane,  then  we  can  neglect  the  trans¬ 
verse  derivatives  in  Eq.  (1)  to  obtain 


£ ( x.y, L, f )  =  £ ( x,y, 0 )exp[ i(t>ix,y,Lj)]  , 
where 


(2) 

_J 


/ dh  f  d^^\E{x,y,0)E(^,r]M^  exp 
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FIG.  3.  Natural  logarithm  of  the  transmission  T  vs  scale 
reading  S  (on  Cl)  for  five  concentrations  of  the  latex  suspension. 
The  particle  diameter  is  1.09  pm  and  the  concentrations  are 
given  in  Table  I.  The  transmission  is  defined  as  the  fraction  of 
power  remaining  in  the  Gaussian  beam  (less  background)  after 
it  traverses  the  latex  suspension,  and  is  calculated  as  explained 
in  the  text.  Solid  lines  are  a  best  fit  according  to  Eq.  (6)  with  the 
scale  reading  L,,  corresponding  to  zero  thickness  being  an  ad¬ 
justable  parameter.  The  values  of  the  parameters  are  given  in 
Table  I.  Curves  F  are  typical  limiting  curves  with  off-center  pa¬ 
rameters  used  to  determine  the  error  bounds  given  in  Table  I; 
the  data  points  are  those  of  curve  A  shifted  up  by  two  units;  the 
curve  F  is  a  plot  of  Eq.  (6)  with  a:^  =0.36  pm  ',  L,.  =  1 1  pm, 
I  .,=  —0.8  pm  and  =0.32  pm“ ',  L,.  =  3  pm,  L()  =  3  pm.  The 
point  P  was  obtained  from  the  typical  OMA  trace  shown  in  Fig. 
2.  The  left-pointing  arrows  indicate  the  points  where  the  beam 
size  corrections  estimated  in  the  text  are  100%.  The  right- 
pointing  arrows  indicate  the  points  where  the  background  in¬ 
tensities  B  are  half  of  the  corresponding  peak  intensities. 


<j>(x,y,L,t)=  f  ^be(x,y,z,t)dz  . 

Inc  •'() 

The  transmitted  fraction  is  then  the  functional  projection 
of  the  field  on  the  exit  plane  on  the  field  that  would  exist 
in  the  absence  of  attenuation.  We  can  express  the  (time 
averaged)  transmission  therefore  by 


f  d~  xE*(x,y,Q)E(x,y,L,l) 

J d^x  E*(x,y,0)E(x,y,Q) 

■) 

where  d'x  =dx  dy  and  all  integrals  extend  to  infinity. 
We  will  assume  that  8e(x,/)  is  a  (zero-mean)  Gaussian 
random  process.  The  numerator  in  Eq.  (3)  reduces  to 


^„-^o\)~Ct]x-i\)]dzd; 


(4) 


where  x,|  =  (0,0,zl,  ;  (0,0,^)  and 
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C(|x-^l)  = 


■( b€{x,y,z,t)b€{^,r],^,t )) 
<8e-) 


is  the  equal  time  correlation  function,  and  G  ={(o/2nc)~(b€~).  We  assume  a  correlation  function  of  the  form 
C(  ix  — 1;|  )  =  exp(  —  jx  — /L;, );  then  the  numerator  becomes 


J d'x  'c|£(x,j>,0)£(5,7;,0)!'exp 


—  G  r^[exp(  —  |z  — ^|/£^. )  — exp(  —  |x-|:| /£,.  Ijtfz 


(5) 


where  £,.,  the  correlation  length,  is  of  the  order  of  10  /xm, 
much  less  than  the  laser  beam  width  m„  of  0.7  mm.  The 
last  term  in  the  exponent  in  Eq.  (5)  gives  a  correction  to 
the  expression  which  is  obtained  without  it,  that  is  of  the 
order  (L,  /(j,,)’  smaller.  Without  this  term  we  get 

r(£)  =  exp[ —a  ,,(£—£,. +£,e  ^  ^' )]  ,  (6) 

where  a^=2GL^  is  the  attenuation  constant  at  large 
thickness. 

We  note  two  points  about  these  approximate  results. 
First,  because  we  assume  that  be  is  a  Gaussian  random 
function  in  space,  the  average  beam  transmission  in  Eq. 
(3),  with  Eq.  (2)  for  the  exit  beam,  implies  that 

((A-)  =  -ln7.  (7) 


Of  course  this  phase  fluctuation  is  only  within  the  exit 
beam  pattern  itself,  Eq.  (2)  does  not  include  the  scattered 
light  at  the  exit  plane.  Second,  when  the  beam  is  much 
wider  than  a  correlation  length,  Eq.  (5)  gives  for  the  rms 
fluctuation  in  Eq.  (7) 


(<b-}=2GL  f' 

•'  (I 


C{z)dz  . 


(8) 


For  L  «£, ,  (S')  approaches  GL',  whereas  for  L  »£,, 
(6')  approaches  a  ^  L.  This  result  follows  from  the  fact 
that  the  phase  d  in  Eq.  12)  undergoes  a  random  walk  as  it 
propagates  inside  the  scattering  medium.  In  particular, 
the  exponential  decay  of  the  transmission  at  large  thick¬ 
ness  is  analogous  to  the  fact  that  the  mean-square  dis¬ 
placement  of  a  Brownian  particle  Is  proportional  to  time 
at  long  times. 


IV.  RESULTS 

The  results  of  fitting  Eq.  (61  to  our  data  are  shown  as 
the  solid  lines  in  Fig.  3.  The  best  parameter  values  ob¬ 
tained  from  this  fit  are  summarized  in  Table  I.  We  see 
that  as  the  particle  density  increases,  the  attenuation  con- 


I - 

slant  a.,,  and  the  correlation  length  increase.  The 
scattering  cross  section  represented  by  the  ratio  a  /N 
does  not  show  any  significant  change  with  density  for 
samples  A,  B,  and  C,  but  decreases  for  the  higher  densi¬ 
ties.  The  error  bounds  given  in  Table  I  are  obtained  by 
simultaneous  variations  of  a  .^  ,  £,.,  and  Lq  from  their  best 
fit  values  until  only  half  of  the  data  points  intercept  the 
theoretical  fit.  An  illustration  of  two  such  limit  curves 
are  shown  by  curves  labeled  F  in  Fig.  3,  and  which  are 
shown  with  the  same  data  as  in  curve  A,  but  shifted  by 
two  units  vertically  for  clarity. 


V.  DISCUSSION 

We  have  assumed  that  the  dielectric  function  8e(x,t)  is 
real.  This  approximation  is  valid  if  the  absorption  of 
beam  energy  by  the  medium  does  not  alter  the  beam 
shape  (we  have  verified  that  it  does  not)  and  if  the  at¬ 
tenuation  arising  from  absorption  is  much  less  than 
or  £,  ',  whichever  is  least.  The  following  experimen¬ 
tal  check  shows  our  nonabsorbing  approximation  to  be 
valid,  even  had  we  used  strong  515  nm  beams  where  ab¬ 
sorption  is  presumably  greater.  We  attempt  to  measure 
the  absorption  constant  of  the  latex  samples  using  the 
thermal  blooming  effect.  We  assume  our  beam  has  a 
transverse  intensity  profile  of  the  form  Hr) 
—  — 2r‘/W|",i,  which  gives  rise  to  an  easily  calcul¬ 

able  transverse  thermal  gradient  in  the  material  due  to  lo¬ 
cal  heating.  The  thermal  gradient  creates  a  refractive  in¬ 
dex  gradient,  thereby  creating  a  thermal  lens  of  focal 
length  F.  Assuming  that  the  beam  length  £  in  the  sample 
is  small  compared  with  both  the  attenuation  length  and 
the  focal  length  of  the  thermal  lens,  and  that  heat  is  only 
being  conducted  in  the  transverse  direction,  Gordon 
cl  al.  showed  that  for  the  above  beam  profile,  this  focal 
length  is  given  by  £  =  77A./((j|V(n„/nF||£//').’  Here: 
k  thermal  conductivity  (J  cm  '  s  'K  '),«  refractive 
index,  absorption  constant  (cm  '),  _  power  of 


TABLE  I.  Mcasurvil  sample  densities  .V  and  best  parameters  for  fitting  Eq.  (6)  to  data  of  Fig.  .E  The 
symbols  are  explained  in  the  text. 
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incident  beam  (If'),  n'  =  [£//j /d(temperature)]  at  room 
temperature  (K  ’),  and  Z,,  and  F  are  in  cm.  Howev¬ 
er,  to  increase  sensitivity,  we  make  the  suspension  that 
we  study  longer  than  the  attenuation  length.  If  the  other 
assumptions  still  hold,  the  expression  of  Gordon  et  al. 
need  only  be  multiplied  by  (lnr)/(7  — 1)  for  this  case, 
provided  that  heating  due  to  the  scattered  light  can  be 
neglected.  Therefore  we  have  been  able  to  set  a  lower 
limit  on  as  follows. 

We  focused  an  argon  laser  (A.  =  515  nm)  to  a  minimum 
beam  cross-sectional  diameter  of  3.6  /am  in  the  suspen¬ 
sion  (resulting  in  a  maximum  intensity  of  1.7  MW/cm’) 
and  were  unable  to  detect  any  thermal  lensing  effect  on 
the  beam.  We  used  a  wavelength  different  from  our  pre¬ 
vious  measurement  so  as  to  have  much  higher  power. 
We  employed  a  sample  like  B  in  Table  I,  for  which  the 
transmission  T  equals  0.125  at  a  thickness  Z.=8  /im. 
With  this  geometry  there  will  be  a  significant  fraction  of 
the  heat  energy  being  transported  along  the  beam  axis. 
This,  we  estimate,  would  not  alter  our  limit  on  by 
more  than  a  factor  of  3.  The  transmitted  beam,  collimat¬ 
ed  by  a  lens  of  focal  length  15  cm  placed  15  cm  from  the 
suspension,  then  passes  through  a  second  identical  lens 
placed  5  cm  from  the  first  lens,  and  impinges  on  the 
OMA  camera  16.5  cm  from  the  second  lens.  There  the 
( 1/fl  beam  diameter  was  observed  to  deviate  less  than 
from  the  low  power  value  (  ^0. 7  mm)  at  all  powers 
up  to  0.7  W.  Using  a  Gaussian  beam  propagation 
analysis  we  obtain  \F^  >0.44  mm  from  this  observation. 
With  this  limit  for  F,  and  the  following  parameter  values: 
Refractive  index  n  =n„;,ui  *-'^57, 

P,,=0.1W,  A'-0.597'Wm  'K  ',  and' 

X  10  ■'  K  '  (the  value  for  water  alone);  we  obtain  the 
satisfying  result  that  the  attenuation  due  to  absorp¬ 
tion  is  less  than  4.5  mm  '. 

We  checked  that  the  transmission  of  the  latex  suspen¬ 
sion  is  independent  of  intensity  so  that  no  nonlinear  opti¬ 
cal  effects  (e.g.,  electrostrictive  effect)**  are  impt)rtant.  We 
also  checked  the  absorption  constant  of  the  dye  solution 
after  the  five  sets  of  measurements,  to  find  that  it  differed 
by  at  most  0.2%.  The  absorption  constant  (0.501  mm  ') 
of  the  dye  used  is  much  smaller  than  the  itienuation  con¬ 
stant  of  the  latex  suspensions  so  that  the  relative  error  in 
the  scale  reading  of  the  dye  solution  is  negligible.  The 
background  intensity  is  much  weaker  than  the  transmit¬ 
ted  intensity  at  transmissions  above  the  right-pointing  ar¬ 
rows  in  Fig.  3.  Therefore,  we  rule  out  multiple  scattering 


as  a  contributor  to  this  nonexponential  behavior.  To 
check  whether  this  behavior  is  a  surface  effect  due  to 
clustering  of  particles  on  the  glass  windows  of  the  cell,  we 
studied  the  sample  under  a  microscope.  We  could  just 
resolve  individual  spheres  at  an  edge  where  they  were 
pinned.  Their  Brownian  motion  was  too  fast  to  render 
them  visible  elsewhere;  no  stationary  spheres  or  clusters 
could  be  seen.  In  fact  we  could  see  no  light  scattering 
centers  at  the  window.  As  another  check  we  performed 
the  same  measurements  with  a  concentrated  dye  solution 
(Oj  ^83  cm  ')  in  cell  Cl,  and  found  perfect  exponential 
decay.  To  test  the  sensitivity  of  our  results  to  the  exact 
form  of  the  correlation  function,  we  performed  a  curve 
fitting  with  another  correlation  function  of  the  form 
C'(  lx  — ^1  )  =  exp(  — 7r|x  — ^|■/4L^’)  which  gives 

7'(L)  =  exp|  — «  ,  [erf(x)L  — 2L,  ( 1  — c  '')/7r]l  ,  (9) 

where  x  =  v  ttL/IL^..  In  our  curve  fitting,  we  made  the 
approximation  erf(.v)  =  tanh(5.v/4),  correct  to  within 
l%f.  A  summary  of  the  best  parameters  and  their  fitting 
limits  is  given  in  Table  II. 

Our  theoretical  treatment  neglected  the  transverse 
derivatives  in  Eq.  (1)  which  govern  not  only  the 
diffraction  of  the  Gaussian  beam,  but  the  diffraction  of 
the  scattered  light  as  well.  We  have  redone  our  analysis 
treating  these  derivatives  in  perturbation  theory.  The  re¬ 
sult  is  a  fractional  correction  of  the  order  of 
(Ao,  to  Eq.  (6)  for  the  transmission  T.  This  cor¬ 

responds  to  a  maximum  correction  of  about  309{-  at  the 
extreme  low  transmission  points. 

We  also  estimated  the  “beam-size”  correction  which, 
we  noted  previously,  arises  because  we  neglected  the  last 
term  in  the  exponent  in  Eq.  (5).  Using  the  expression 
above  Eq.  (9)  as  the  correlation  function,  we  found  the 
lowest-order  fractional  correction  to  be  of  the  order 
L,"/(wornn7'i )  which,  for  our  geometry  is  much  larger 
than  the  diffraction  corrections.  In  Fig.  3,  the  points 
where  this  beam  size  correction  is  100%  is  indicated  by 
left-pointing  arrows. 

At  small  thickness  in  sample  B,  the  scattered  light  is 
:-.trong]y  peaked  within  a  3°  radius  cone  about  the  for¬ 
ward  direction  where,  presumably,  it  consists  mainly  of 
singly-scattered  photons.  This  structure  factor  [which  is 
approximately  proportional  to  the  scattered  intensity 
liO)]  in  the  near  forward  direction  is  well  approximated 
by  a  Gatissian  function  plus  a  constant  and  the  correla- 


TABl.F  II.  Measured  sample  densities  ,V  and  best  parameters  fur  litting  Eq.  (9l  to  data  of  Fig.  .F 
I  he  symbols  are  explained  m  the  text.  Measurements  of  the  structure  factor  in  the  near  forward  direc¬ 
tion  gives  a  correlation  length  of  X  2  /on  for  sample  /i  and  7  2  //m  for  a  sample  with  particle  dcnsitv 
one-lilih  of  that  of  sample  H 
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tion  function  obtained  by  Fourier  transform  suggests  [ac¬ 
cording  to  the  correlation  function  above  Eq.  (9)]  a  corre¬ 
lation  length  in  the  scattering  medium  of  8±2  fim,  in 
agreement  with  the  we  find  from  our  model. 

There  are  three  types  of  forces  contributing  to  the  in¬ 
teraction  between  the  particles.  The  first  is  the  hard 
sphere  collision  between  the  paiiicles  wmch  gives  a 
correlation  length,  according  to  the  Percus-Yevick  ap¬ 
proximation,*^  which  is  approximately  that  of  the  particle 
diameter.  This  is  different  from  what  we  observe  in  the 
near  forward  direction.  The  second  is  the  hydrodynamic 
interaction  which  is  dependent  on  the  density  of  the  par¬ 
ticles.  We  have  measured  the  structure  factor  in  the  near 
forward  direction  for  a  suspension  with  particle  density 
one-fifth  of  that  of  sample  B  and  obtained  the  correlation 
length  as  7±2  /im.  This  suggests  that  the  third  kind  of 
interaction,  the  electrostatic  force,  is  dominant  in  our 
suspensions.  This  also  agrees  with  our  dc  conductivity 
measurements  since  the  measured  conductivity  is  two  or¬ 
ders  of  magnitude  higher  than  that  of  water. 


VI.  SUMMARY 

In  summary,  we  have  observed  experimentally  a  strong 
deviation  from  Bouguer’s  law  of  exponential  beam  at¬ 
tenuation  in  the  low-intensity  limit.  These  data  are  very 
well  fitted  by  a  model  which  assumes  Gaussian  random 
fluctuations  in  the  dielectric  function  that  liave  a  correla¬ 
tion  length  longer  than  the  attenuation  length.  This 
correlation  length  agrees  with  that  estimated  from  the 
angular  spectrum  of  the  scattered  light.  We  verified  that 
multiple  scattering,  thermal  blooming,  surface  clustering 
effects,  and  nonlinear  optical  effects  do  not  contribute  to 
this  deviation  from  exponential  attenuation. 
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Direct  determination  of  eiectron  mobiiity  in  photorefractive  Bii2Si02o 
by  a  hoiographic  time-of-fiight  technique 
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The  development  of  a  photorefractive  grating  created  by  intersecting  30  ps  (532  nm) 
beams  in  a  well-characterized  Bii2SiO20  crystal  (in  a  static  field  around  1  kV/cm) 
unambiguously  reveals  the  mobility  of  photoexcited  electrons  to  be  0.24  ±0.07  cm^/V  s 
through  what  is  essentially  a  “time-of-flight”  measurement. 


Electrical  measurements  of  the  mobility  of  photoex¬ 
cited  carriers  in  photorefractive  insulating  crystals  have 
been  inconsistent  and  difficult  to  interpret,  in  part  because 
of  electrode  potentials,  surface  currents,  and  hot  carriers. 
Mobility  values  from  5  X 10  ~  ’  to  3  cm^/V  s  have  been 
reported  for  undoped  n-type  cubic  Bii2SiO20  (n-BSO) 
from  the  photocurrent  measurements. Here  we  report 
holographic  time-of-flight  measurement  of  electron  drift 
mobility  fi,  finding  0.24  ±0.07  cm^/V  s.  The  only  previous 
holographic  time-of-flight  measurement  of  the  mobility  of 
photorefractive  carriers  known  to  us  is  that  of  Pauliat  et 
al.*  who  found  values  around  10  “  ^  cm^/V  s  for  the  carri¬ 
ers  of  unknown  charge  sign  in  three  cubic,  iron-doped 
Bi|2GeO20  crystals.  The  technique  we  report  here  is  similar 
to  that  employed  in  Ref.  4  except  that  our  carriers  drift 
under  the  influence  of  a  static  applied  voltage  Vq  rather 
than  an  alternating  one,  and  we  excite  the  carriers  with  an 
essentially  instantaneous  532  nm  pulse  (ss30  ps,  1  fjJ/ 
cm^),  rather  than  with  a  cw  beam. 

Described  briefly,  our  technique  for  observing  drift 
mobility  works  as  follows.  A  density  “grating” 
n(x,r)(a£«o(f) -i- Re[/i|(0exp(/far)])  cm“^  of  photoex¬ 
cited  carriers  (electrons  here)  is  initially  superposed  on  the 
grating  of  ionized  deep  traps  from  which  the  carriers  were 
instantaneously  excited.  This  grating  n{x,t)  drifts  subse¬ 
quently  under  the  influence  of  the  xr-directed  static  uniform 
electric  field  Eq.  Before  significant  diffusion  takes  place,  we 
observe  the  rise  and  fall  of  the  amplitude  Ei  of  the  Cou¬ 
lomb  electric  field  grating  (  Re[E|(r)exp(z7cx)])  caused  by 
the  movement  of  the  electron  grating  from  coincidence,  to 
anticoincidence,  and  back  to  coincidence  with  the  station¬ 
ary  grating  of  ionized  traps.  (£|  <£0)  The  time  for  this  rise 
of  the  Coulomb  field  grating  (observed  here  to  be 
10“ ’ s)  is  evidently  ir/{  fiEffc),  where  the  grating  wave 
vector  k  equals  [2ir/grating  period].  Therefore,  we  are  jus¬ 
tified  in  calling  our  measurement  a  direct  "time-of-flight” 
mobility  measurement.  We  monitor  £'|(/)  by  observing  the 
Bragg  scattering  of  a  weak  laser  beam  from  the  refractive 
index  grating  An(x,r)  that  £|  causes,  via  the  electro-optic 
effect:  An(x,t)  =  Re[  —  (l/2)nV£|  exp(/Ayt)].  Here  n  is 
the  linear  refractive  index  and  r  is  the  effective  electro-optic 
coefficient  for  the  geometry.  We  also  monitor  Eq  via  the 
electro-optic  effect.  In  the  following  sections,  we  first  out- 
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line  a  simple  theory  which  includes  the  effects  of  electron 
diffusion  and  recombination,  and  also  the  effects  of  the 
spatial  variation  of  the  mobility.  We  then  review  the  mo¬ 
bility  experiments  which  we  performed  as  described  above 
on  a  5x5x5  mm^  single  crystal  of  nominally  undoped 
n-type  Bii2SiO20  of  high  optical  quality,  grown  by  Sumit¬ 
omo  Inc.,  and  well  characterized  by  the  various  experi¬ 
ments  described  in  Ref.  5  where  this  crystal  was  referred  to 
asSUl. 

We  may  obtain  a  more  precise  expression  for  the  form 
of  !£](/)  1^,  which  the  Bragg-scattered  633  nm  beam  mon¬ 
itors,  by  using  the  Poisson  equation 

ikEi=eiNi-ni)/e,  (I) 

where  e  is  the  dielectric  constant  of  the  crystal,  e  is  the 
(positive)  electronic  charge,  and  Ni(t)  cm  "  ^  is  the  (com¬ 
plex)  amplitude  of  the  spatially  sinusoidal  density  grating 
of  ionized  traps.  Since  we  excite  many  fewer  electrons 
( ~  10’^  cm  “  ^)  than  the  dark  density  of  ionized  traps 
( =£  10'*  cm  “  ^  from  Ref  5)  we  may  assume  that  A^i  obeys 
a  simple  equation  for  direct  recombination: 

-  n,/T,  (2) 

where  r  is  the  usual  electron  recombination  time,  which  we 
have  measured  by  direct  observation  of  photocurrent  decay 
to  be  80±5  fis.  Initially  Ai(0)  =n|(0)  (before  the  elec¬ 
trons  have  moved)  and  £j(0)  =0.  We  assume  that  the 
subsequent  electron  current  density  j(x,t)  is  in  the  x  direc¬ 
tion  and  is  related  to  the  electron  density  n  by  the  usual 
equation  representing  drift  and  diffusion: 

j{x,t)  =fien[Eo  -(-  Re(£,e''“)  ]  -f  (ikgTdn/dx.  (3) 

Here  is  the  Boltzmann  constant  and  T  the  temperature. 
Our  experiments  are  in  the  limit  of  low  photoexcitation 
where  we  see  only  the  terms  in  j  which  are  linear  in  the 
energy  flux  U  J/cm^  of  the  532  nm  excitation  pulse.  Since 
both  £)  and  n  are  linear  in  U  in  this  limit,  we  do  not 
observe  the  second  term  of  (3).  The  unimportance  of  this 
term  is  confirmed  by  our  theory  which  predicts  |  £0 1  <  |  £]  | 
for  our  experimental  conditions:  £/=:10~*  J/cm^.  Omit¬ 
ting  this  small  term  allows  us  to  eliminate  j  from  (3)  by 
using  the  continuity  equation 

dj/dx=e(dn/dt n/r)  (4) 

to  obtain  for  the  complex  amplitude  n\  of  the  density  grat¬ 
ing 
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n,(f)=/i,o^  (5) 

where  «,o  =  «i(0)  and 

r  =  r'  +  '  —  ik[x(Eo  +  ikkgT/e).  (6) 

With  this  in  (2)  we  have  immediately  for  the  desired  Cou¬ 
lomb  field  grating  amplitude 

(7) 

where  =  —  nen]o(Eo  -|-  ikkgT/e)/i€r). 

In  our  experiment  we  arranged  the  grating  period  k  to 
be  670  cm  “  Then  kkgT/e  is  17  V/cm  which  is  much 
smaller  than  the  fields  Eq  (300-2500  V/cm)  which  we 
applied  to  the  crystal.  Also  {k^Lr)  ~ '  is  77  V/cm  in  our 
experiment.  Therefore,  T  is  nearly  (  —  ikfiEg)  which  cre¬ 
ates  the  sinusoidal  rise  and  fall  of  l^ij^  explained  in  our 
introduction  above.  The  relative  dielectric  constant  e/tg  of 
our  BSO  sample  (called  SUl  in  Ref.  5)  was  56  so  that,  for 
the  approximately  1.6xl0’^  cm“^  electrons  initially  ex¬ 
cited  by  the  532  nm  intensity  grating,  we  expect  from  (7) 
that  I  £,  I  will  have  a  maximum  value  near  80  V/cm  and 
will  not  depend  much  on  any  variation  in  Eq,  or  r  that 
may  occur  in  our  experiment. 

We  have  actually  observed  more  damping  in  the  oscil¬ 
latory  behavior  of  |  |  versus  time  than  is  predicted  by 

(7).  We  can,  however,  obtain  an  excellent  agreement  with 
experiment  which  we  document  below  if  we  assume  that 
the  mobility  n  varies  randomly  with  position  x  in  the  crys¬ 
tal  on  a  scale  larger  than  a  grating  spacing  but  with  a 
correlation  length  much  smaller  than  any  dimension  of  the 
Bragg-scattering  region.  The  total  Bragg-scattered  (phase- 
matched)  633  nm  signal  is  essentially  a  coherent  superpo¬ 
sition  of  the  radiation  generated  by  the  electro-optic  effect, 
and  hence  has  a  power  proportional  to  \S  dx^ 
Ep{x)E^(x)  1^,  where  Ep  is  the  amplitude  of  the  633  nm 
beam  scattered,  and  the  integral  is  taken  over  the  volume 
of  the  crystal.  We  assume  that  filx)  varies  through  the 
variation  of  ^i{x)  and  T  in  (7),  and  then  average  over  the 
deviations  of  ^  from  its  average  Ji.  For  simplicity  we 
assume  that  6/r  is  a  Gaussian  random  variable  and  so  has 
the  property  that  (exp(/)36/i))  =  exp(  —  ( l/2)/3’<6^r) ). 
Neglecting  the  variations  in  E^  arising  from  5^i  (for  rea¬ 
sons  explained  above)  we  evaluate  the  foregoing  space  in¬ 
tegral  using  (7)  for  £,  to  find  directly  that  the  intensity  Ig 
of  the  Bragg-scattered  signal  is  predicted  to  obey 

-exp[  -  (PJl  t)~C^/2  +  iPJI  t  ~  t/T]\^, 

(8) 

where  P  =  kE^  +  ik^kgT/e  and  =  {(d^i)') /Ji'  is  a  pa¬ 
rameter  describing  the  fractional  fluctuation  of  mobility 
throughout  the  crystal.  We  find  below  that  taking  C  to  be 
0.42  makes  (8)  to  fit  all  our  data  to  within  experimental 
error.  This  is  not  an  unreasonable  amount  of  fluctuation  in 
mobility,  given  the  uncontrolled  nature  of  the  electron  do¬ 
nors  and  of  the  impurities  in  the  growth  process.  The 
infinite-time  asymptote  I ^  is  proportional  to  \E^  |'. 

Our  experimental  arrangement  is  diagrammed  in  Fig. 
1.  The  beam  expander  B-E  and  beamsplitter  BSI  form  two 


FIG.  I.  Experimental  arrangement  for  the  holographic  time-of-fught 
measurement. 

beams  of  30  ps  532  nm  pulses  labeled  S]  and  Si  which  are 
shuttered  by  Sh|  and  Shj.  The  beam  Sj  excites  the  initial 
grating  of  photoelectrons,  its  flux  having  been  spatially 
modulated  to  equal  IJ[\  +  cos(A:x)]  in  the  BSO  sample,  by 
transmission  through  the  Ronchi  ruling  R  and  the  spatial 
filter  (compromising  lenses  L\  and  L2  separated  by  aper¬ 
ture  A ) .  The  0  and  1  orders  were  used  in  our  experiment  to 
make  k  =  670  cm"  '  (a  grating  period  of  93  /um).  The 
experimentally  employed  range  of  U  was  between  0. 1  and 
5  /iJ/cm^.  The  second  532  nm  beam  Si  was  used  to  erase 
previously  written  photorefractive  gratings  before  the  ex¬ 
citation.  Both  beams  had  an  intensity  that  was  uniform 
over  the  He-Ne  probe  beam  to  within  5%. 

The  variation  in  time  of  the  photorefractive  gratings  in 
the  BSO  sample  was  monitored  by  the  (1.2  mm  diameter) 
633  nm  He-Ne  laser  beam  controlled  by  an  acousto-optic 
modulator  ( AOM),  and  linearly  polarized  by  £,  at  such  an 
angle  as  to  make  its  Bragg-diffracted  signal  orthogonally 
polarized  to  the  nearby  directly  transmitted  beam.  The  po¬ 
larizer  Pi  was  then  effective  in  eliminating  ordinary  scat¬ 
tered  light  from  the  Bragg-diffracted  signal  recorded  by  the 
photomultiplier  (PMT). 

The  voltage  which  was  applied  to  the  BSO  sample 
to  create  the  drift  in  the  x  direction  of  the  photoexcited 
electrons,  was  carefully  coordinated  in  time  with  the  shut¬ 
ters  Sh,  and  Sh2  and  the  acousto-optic  modulator  to  min¬ 
imize  the  residual  space  charge.  The  polarization  of  the 
primary  633  nm  beam  that  was  transmitted  through  the 
sample  was  measured  in  order  to  determine  the  applied 
electric  field  f,,  via  the  electro-optic  effect  as  described  in 
Ref.  6.  We  found  that  Eq  was  spatially  uniform  and  typi¬ 
cally  70%  of  Viy/d,  where  d  was  the  distance  between  (5 
mm)  the  coated  electrodes.  We  believe  that  the  accuracy  of 
the  electric  field  measurement  was  ±25%.  The  direction 
(x)  of  Eq  was  parallel  to  the  ( 1.1,0)  crystal  axis.  The  532 
nm  beams  were  parallel  to  the  ( l.T.O)  axis. 

The  diffracted  633  nm  intensity  versus  time  averaged 
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FIG.  2.  Bragg-scattered  intensity  of  the  633  nm  beam  vs  time  /  from  the 
BSO  sample  after  the  picosecond  photoexcitation  pulse  when  four  dif¬ 
ferent  voltages  were  applied  across  the  sample.  Also  two  best-fit  curves 
are  shown. 

over  about  50  shots  is  shown  in  Fig.  2  for  four  different 
applied  voltages  to  the  sample  “SU 1  ”  of  n-type  BSO.*  The 
two  solid  curves  drawn  over  each  data  set  correspond  to 
two  choices  for  the  set  of  parameters  (C,\PJi\,I^)  in  a 
plot  of  (8).  The  independently  observed  recombination 
time  r  =  80  /iS  and  a  small  parameter  d  in  PJI  =  \ 
(deduced  from  the  electric  field  measurements)  were  used 
in  all  of  the  solid  curve  plots.  The  solid  curve  with  the 
lowest  maximum  corresponds  to  C  =  0.42  in  each  case 
along  with  the  least-squares  values  of  |  Pji\  and  listed 
in  Table  I.  The  other  solid  curve  is  for  C  =  0.35  in  each 
case,  value  which  increases  the  rms  deviation  of  the  data  by 
20%.  The  errors  listed  in  Table  I  for  (  PJI)  and  reflect 
the  new  best  values  when  C  =  0.35  is  used.  The  variation  in 


is  mainly  caused  by  the  differences  in  the  excitation 
energy  fluxes  (/^  proportional  to  |  C/l"). 

One  sees  that  none  of  the  four  fields  gives  a  value  for 
the  mobility  Ji  that  deviates  by  more  than  5%  from  their 
mean  value  0.24  cm’/V  s.  We  feel  that  the  uncertainty  in 
the  true  value  for  this  sample  arises  mainly  from  our  esti¬ 
mated  25%  uncertainty  in  the  absolute  value  of  Eq  which 
we  calculated  from  the  polarization  shift  of  the  633  nm 
beam  (caused  by  the  electro-optic  effect)  assuming  the 
electro-optic  coefficient  of  our  sample  to  be  the  gener¬ 
ally  accepted  value  4.4  pm/V  for  BSO  at  633  nm.’  There¬ 
fore  we  have  assigned  an  estimated  uncertainty  of  ±0.07 
cmVv  s  to  our  mobility  determination. 

In  any  case,  we  note  that,  had  we  simply  measured  the 
time  for  the  signal  to  reach  its  peak,  and  equated  that 
time  to  n/iJIEQk)  we  would  have  obtained  the  same  value 
for  ju  to  within  5%  in  each  case.  Our  more  careful  analysis 
of  the  data  should  not  cloud  the  fact  that  we  have  essen¬ 
tially  a  “time-of-flight”  measurement  of  the  mobility  here. 
The  most  interesting  insight  to  come  from  our  data  anal¬ 
ysis  is  that  our  data  fit  the  theoretical  expression  (8)  much 
better  if  we  assume  that  the  macroscopic  mobility-field 
product  (/iEq)  fluctuates  by  ±40%  from  place  to  place  in 
the  crystal.  By  “macroscopic”  we  mean  to  imply  that  fluc¬ 
tuations  on  a  scale  smaller  than  the  grating  period  do  not 
affect  our  results.  Since  we  could  not  see  variations  in  Eq 
larger  than  ±5%  using  the  electro-optic  effect,  we  feel 
justified  in  assuming  that  the  mobility  itself  was  varying 
from  place  to  place  in  the  crystal  by  roughly  ±40%. 

We  note  also  that  the  value  for  the  product  (fir)  for 
the  electrons  in  SU  1  was  inferred  from  the  k  dependence  of 
grating  erasure  rates  in  Ref.  6.  With  the  value  t=  80±5 
fis  which  we  have  measured  by  photocurrent  decay,*  we 
would  predict  /x  =  0.1 7 ±0.1  cm’A' s,  a  value  consistent 
with  our  more  accurate  finding  here.  If  we  assume  the 
quantum  efficiency  for  photoelectrons  to  be  =85%  as  es¬ 
timated  in  Ref  6,  then  our  new  mobility  value  correctly 
predicts  the  absolute  value  of  current  in  the  photocurrent 
experiment  with  a  30  ps  532  nm  light  pulse.*  The  variation 
of  in  Table  I  is  also  consistent  with  that  expected  using 
the  parameters  of  Ref.  6  in  (8). 

This  research  was  supported  by  the  U.  S.  Air  Force 
Office  of  Scientific  Research  under  contract  No.  F  49620- 
88-C-0027  and  by  the  National  Science  Foundation  under 
grant  No.  ECS-8821507. 


TABLE  I.  Spatial  average  JI  of  photoelectron  mobility  as  obtained  by 
adjusting  \pfl\  and  /,  to  obtain  the  best  fit  of  (8)  to  the  data  of  Fig.  2. 
for  four  values  of  applied  voltage  V,,  to  Bi,;SiO;o  sample  SUI.  The  least- 
squares  error  value  C  =  0.42  was  used  in  (8)  for  the  computation. 


E, 

(kV) 

V 

(pj/cm’) 

F 

(kV/cm) 

o - 

/, 

(10 

Jt 

(cm’/Vs) 

0.25 

0.78 

0.31 

50.8  ±0.5 

80±2 

0.245 

0,50 

0.73 

0.67 

105  ±2 

90  ±3 

0.2.34 

1  00 

093 

1.38 

2.30  ±5 

157  ±7 

0.249 

1  50 

099 

2.17 

362  ±  12 

183  ±9 

0.249 

■“These  values,  measured  as  explained  in  text,  have  relative  uncertainly 
5%  and  absolute  uncertainly  of  25'5r. 
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RESUME  Nous  prtsentons  une  tnesure  holographique  du  temps  de  vol  des  poneuis 
de  charge  dans  les  mat6riaux  photortfractifs  qui  permet  d'acc6der  directement  k  leur 
mobility.  Nous  I'avons  exp6riment6e  sur  un  echantillon  de  Bii2SiO20  pr6alablement 
caracterisd  par  diH'ifrentes  techniques.  Nos  rdsultats  montrent  que  la  vitesse  d* 
cntrainement  des  poneurs  est  proportionnelle  au  champ  electrique  qui  leur  est  appliqu^ 
(de  2(X)  k  20(X)  Vcm'^  environ).  On  trouve  alors,  ind6pendamment  de  tout  autre 
paramitre  du  matdriau,  une  mobility  de  0.24±0,07  cm^V'^s'*.  Nous  prisentons  aussi 
des  mesures  de  photocourants  transitoires  confiimant  ce  lesultat. 

ABSTRACT  We  present  a  holographic  time-of-flight  techniq  le  to  access  directly  the 
mobility  of  the  photoexcited  charge  carriers  in  photorefractive  materials.  We 
experimented  it  on  a  previously  well  characterized  sample  of  Bi)2SiO20-  Our  data 
show  a  linear  dependence  of  die  drift  velocity  versus  applied  electric  Held  in  our 
experimental  range  (-200  to  2000  Vcm'*).  We  find,  independantly  of  any  other 
material  parameter,  a  mobility  of  0.24±0.07  cm^V'^s'l  which  is  consistent  with  other 
transient  photocunent  experiments  also  presented  here. 


The  photorefractive  effect,  known  to  allow  near  real  time  image  processing  and  phase 
conjugation  at  veiy  low  illumination  has  been  intensively  investigated  in  semi-insulating 
materials  like  BaTiOs  ,  LiNbOs,  Bii2SiO20  (BSO).  There  is  now  a  general  agreement  on  the 
mechanism  of  the  effect:  The  light  induced  generation  of  free  carriers  from  photorefractive  sites, 
followed  by  their  transport  and  recombination  leads  to  a  charge  redistribution  and  a  space  charge 
Coulomb  field  which  is  turned  into  an  index  variation  through  the  linear  electrooptic  effect. 
However,  neither  the  nature  of  the  photorefractive  sites  nor  the  charge  transport  mechanism  are 
yet  fully  understood.  The  study,  mostly  in  BSO,  of  the  transient  photocurrents  generated  by 
laser  pulses  under  an  applied  dc  electric  field  provides  estimates  for  the  mobility  of  the  charge 
carriers  ranging  from  5  10*5  to  3  cm2/(Vs)*'2.  Values  have  also  been  obtained  for  the  diffusion 
length  (which  caracterizes  the  product  of  the  mobility  4  by  the  lifetime  X  of  the  charge  carrier ), 
from  beam  coupling  or  decay  of  the  diffraction  efficiency  of  a  photoreftactrivc  grating  under 
cw-illumination^.  Using  those  values,  an  estimate  of  the  mobility  in  BSO,  as  high  as  SO 
cm2/(Vs),  was  obtained  from  the  build-up  of  a  grating  of  short  period  after  a  28ps  pulse^- 
Recently,  Pauliat  et  al.  obtained  a  value  4  10'3  cm2/(Vs)  for  a  Bii2GeO20  (BGO)  cryst^  from 
the  time  dependence  of  the  beam  coupling  gain  in  a  two  wave  mixing  experiment  with  a 
rectangular  A.C.  field  applied  to  thecrystal^.  To  our  knowledge,  this  was  the  first  published 
time-of-flight  technique.  These  different  methods  give  very  low  values  for  the  mobility.  Using 
them,  the  mean  free  path  of  a  charge  carrier,  estimated  in  a  simple  model  where  its  movement  is 
limited  by  collisions  would  be  a  fraction  of  an  angstrom  (much  smaller  than  any  interatomic 
distance  in  the  crystal).  The  technique  we  describe  here  is  similar  to  that  published  by  Pauliat 
except  that  our  carriers  drift  under  the  influence  of  a  static  electric  field  after  having  been 
photoexcited  by  an  essentially  instantaneous  532  nm  pulse  (  *  30  ps,  IpJ/cm^)  rather  than  by  a 
cw  illumination. The  success  of  this  versatile  technique  requires  the  possibility  for  the 
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photoexcited  carrier  to  drift  at  least  one  holographic  grating  period  before  recombining.  The 
pulse  energy  should  be  able  to  induce  a  Coulomb  field  grating  strong  enough  to  be  observed 
through  the  electrooptic  effect  but  weak  enough  to  avoid  any  saturation  of  ^e  excitation  and 
transport  presses.  From  a  standard  photocurrent  experiment,  we  find  the  experimental 
conditions  in  which  these  requirements  are  fulfilled.  We  report  our  result  in  the  case  of  a 
previously  well  characterized  crystal  and  show  that  our  values  for  the  mobility  are  consistent 
with  the  values  projected  from  our  photoconductivity  data. 

In  section  n,  we  establish  the  principle  of  the  holographic  time-of-flight  method,  within 
the  simplest  band  conduction  model,  in  a  regime  of  low  illumination.  In  Section  m,  we  outline 
the  results  of  a  transient  photocurrent  experiment  that  set  the  proper  experimental  conditions  fcx' 
the  time-of-flight  experiment.  The  apparatus  and  results  of  the  time-of-flight  experiment  are 
described  in  section  IV.  In  section  V,  we  compare  these  results  to  those  implied  from  our 
photocurrent  trteasurements  and  summarize  our  conclusions. 

II-THE  HOLOGRAPHIC  TIME-OF-FLIGHT  METHOD. 

Consider  a  phorotrefractive  crystal,  illuminated  by  the  intensity  pattern 

I  (r.z)=/o(0(l  +  ‘^*{wexp(tl:.r )})  (1) 

resuldng  from  the  interference  of  two  short  laser  pulses  of  duration  tp  .  In  the  simplest  model 
for  the  photorefracdve  effect,  including  a  single  type  of  photorefractive  sites  of  density  No,  and 
a  single  type  of  charge  carriers  (e.g.,  electrons),  we  study  the  dynamics  of  the  ionized  site 
density  Nd‘  ,  the  density  n  of  conduction  band  electrons  and  the  electric  field  E  inside  the 
crystal. 

The  pulse  ouration  Zp  is  assumed  to  be  shorter  than  any  other  caracteristtc  time,  so  that 
neither  recombination  nor  noticeable  charge  movement  occurs  during  Tp:  no  space  charge  field 

builds-up  during  the  pulse.  If  the  energy  of  the  pulse  is  small  enough  {5/Tp«l),  the  initial 
conditions  may  be  wrinen  as  ; 

^  A  -A'a  )(l  +  ^»{mexpikz}) 

(3) 

n(z,Tp}=S/oTp(N  o  expikzj) 


After  the  pulse  (i.c.,  in  the  dark)  the  evolution  of  the  system,  is  controlcd  by  four 
equations.  The  rate  equation  describes  the  recombination  process  : 

where  yis  the  recombination  rate  constant.  Assuming  free  electron  transport  in  the  conduction 
band,  the  current  equation  includes  both  diffusion  and  drift  components  near  thermal 
equilibrium  at  temperature  T, 

j(z,l)  =  n(z,i)e  ^  £(z,t)  +  fikBT  Vn(z,t)  W 

where  e  is  the  charge  (positive)  of  the  electron,  /t  its  mobility  (positive)  and  kg  the  Boltzmann’s 
constant  The  continuity  equation  for  the  photo-excited  electrons  is, 

^—T^—+yNQ{z,t)n{z,t)=-V.j{z,t).  (5) 

at  e 

The  Poisson  equation  describes  the  electric  field  resulting  from  both  electron  and  ionized  donor 
densities 

V£(z,t)=-^(N  p{z,i)-N  ^  -n{z,t))  (6) 

where  e  is  the  permeability  of  the  crystal  and  Na  the  density  of  acceptors  required  by  electrical 
neutrality.  Assuming  a  small  modulation  m,  ^s.  (3)  -  (6)  may  be  linearized  by  writing  the 
spatially  variable  quantities  as  the  sum  of  a  constant  term  and  of  a  small  sinusoidal  one  : 


(7) 


N  q(2J)=N  o(t)  +M»  [N  ,(f  )cx’p{ikz )} 
n  (2  ,r )  =/io(0  {«iU )  txp(ikz )} 

£  (r  ,f )  =  £  0  +  f  1  (^ )  exp(  ikz)} 


where  Eq  is  a  dc  electric  field  applied  to  the  crystal.  A  system  of  linear  differential  equations  is 
obtained.  Here,  we  solve  it,  in  the  case  of  low  energy  pulses  i.e  when  the  conditions 

nQ{Xp)«N  ^  and  /io(/)«n,  =  —  ,  are  satisfied.  Both  no  and  £/  are  then 

/te  e 

small  and  the  terms  in  no  may  then  be  neglected  in  the  equations  for  the  spatially  varying 
components  of  (7).  Two  simple  differential  equations  are  obtained  : 


f 

y^A 

V 

1+^ 

K 

\  J 

-ikiiE  0 

) 

n,(t)  =  0 

3  t 


^  ^  =  yN  ^  -^(i£  +1^  )n,(/) 

a  t  cK^ 


(8) 
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with ; 

They  are  easily  solved  as 

£,(r)=£4l-e''’) 


KgT 


r s r’+t r ■ = yv ^  —  iknE o 

e  r 


(9) 


(10) 


The  electric  field  shows  damped  oscillations  caused  by  the  superposition  of  the 
stationary  ion  grating  upon  a  grating  of  free  electrons  drifting  with  constant  velocity  under  the 
influence  of  the  electric  field  Eq.  The  oscillations  are  associated  to  the  successive  coincidences 
and  anti-coincidences  of  these  two  sinusoidal  distributions  .  They  are  decaying  progressively  by 
diffusion  and  recombination.  If  the  dc  electric  field  applied  to  the  crystal  is  strong  enough  (ks  > 

k  +  i^/k ),  the  imaginary  part  F"  of  the  damping  constant  may  be  smaller  than  the  real  pan  P. 
The  free  carrier  grating  may  "fly"  one  or  several  holographic  grating  periods  before  vanishing. 
.P'  does  not  depend  on  any  other  material  parameters  than  the  mobility  itself.  It  can  be  accessed 


by  diffracting  a  laser  beam  from  the  crystal.  The  diffraction  efficiency  t)  is  proportional  to  the 
square  of  the  amplitude  of  the  Coulomb  grating  £/: 


(11) 


where  qoo  is  the  value  of  the  diffraction  efficiency  after  damping  of  the  oscillations.  This  method 
is  a  direct  illustration  of  the  concept  of  mobility  that  states  that  free  carriers  move  with  a  speed 
proponional  to  the  electric  field.  Measuring  F"  allows  undisputed  determination  of  that  mobility. 
Only,  the  value  of  the  electric  field  inside  the  crystal  needs  to  be  correctly  evaluated.  Before 
presenting  our  experimental  set-up  and  fitting  our  results  to  equation  (1 1),  we  summarize  our 

transient  photocurrent  experiments  that  set  the  proper  conditions  to  obtain  F"  »P. 


ni  TRANSIENT  PHOTQCURRENT  EXPERIMENTS. 

We  studied  the  time  evolution  of  the  photocurrent  induced  by  a  picosecond  laser  pulse  in 
a  5  X  S  X  S  mm^  crystal  by  Sumitomo  named  SUl  in  Ref.3  and  previously  characteriied  in  that 
paper  by  fitting  two-beam  coupling  and  grating  erasure  data  to  the  electron-hole  competition 

model.  These  experiments  determine  the  mobility  x  lifetime  product  px  of  the  carriers  but  cannot 
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access  them  separately.  Our  photocurrent  experiment  will  provide  an  estimate  for  x  and  thereby 
allow  us  to  estimate  the  proper  grating  period  for  the  time-of-flight  measurement 

The  photocunent  i(t)  induced  by  a  picosecond  pulse  was  studied  with  an  electric  field 
(=lkVcm'*)  applied  to  the  crystal.  The  arrangement  is  shown  in  Fig.l.  The  light  source  is  a 
mode  locked  Nd:YAG  laser  from  which  pulses  are  extracted  with  a  repetition  rate  of  5  Hz  using 
a  Pockels  cell.  After  amplification  and  frequency  doubling,  30  ps  pulses  at  532nm  are  obtained. 
An  electromechanical  shutter  is  used  to  select  a  single  pulse  and  the  laser  beam  is  expaixled  for 
uniform  illumination  of  the  whole  crystal.  Two  photodiodes  PDi  and  PDi  monitor  the  incident 
and  transrruned  energies.  Voltages  Vp  up  to  IkV  can  be  applied  to  the  crystal  by  means  of  silver 
paint  electrodes.  The  photocurrent  is  then  measured  using  a  carefully  designed  ciscuit,  with 

short  connections  to  minimize  stray  inductance,  a  large  high  voltage  capacitor  C  and  50£2  coaxial 
cable.  It  has  a  rise  time  of  1  ns  and  less  than  10%  ringing.  When  a  high  resolution  is  needed,  the 
50fl  cable  itself  plays  the  role  of  the  load  resistor  and  is  directly  connected  to  the  50£2  input  of  a 
400MHz  oscilloscope.  For  longer  time  scales  and  high  sensitivities,  we  use  a  lOOkQ  load 
resistor  connected  to  the  IMQinput  of  the  same  oscilloscope. 


Fig.l  Experimental  set-up  for  the  measurement  of  transient  photocurrents  in  photorefractive 
crystals. 


0  20  40  60  80  0  1 00  200  300  400 
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Fig.2  Transient  photocurrent  in  crystal  SUI  after  a  single  30  ps  laser  pulse  a)  fast  exponential 
decay  ti=26±3  ns  -  b)  Slower  exponential  decay  t2=80i5ps  oberved  with  lower  resolution  set¬ 
up  ( see  text). 
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At  low  fluences  (S  SO^iJ/ctn^)  the  photocurrent  shows  essentially  two  exponential 

decays  having  rime  constants  of  ti=26±3  ns  and  t2=80±5^s  (Fig.2). 

Although  the  photocurrent  fs  much  smaller  during  the  larger  rime  scale,  about  100  rimes 
more  charges  flow  inside  the  crystal  during  that  period  of  rime  than  during  the  initial  current 
pulse.  The  photocurrent  data  can  be  conveniently  expressed  in  the  form  of  the  average 
magnitude  of  the  distance  xj,t)  travelled  by  a  charge  carrier : 

(12) 

nsl 

Where  N  is  the  number  density  of  free  charge  carriers  at  the  end  of  the  laser  pulse,  d  is  the 
distance  between  the  electrodes  and  XnU)  the  distance  travelled  by  a  charge  in  the  direction  of 
the  applied  electric  field  from  its  original  position,  up  to  time  t. 

We  computed  Xa(0  by  integrating  the  photocurrent.  N  was  estimated  assuming  a  unit 
quantum  efflciency  as  in  Ref.3,  taking  into  account  20%  reflections  at  the  surfaces  of  the  crystal 
and  using  an  absorption  coefficient  of  0.74±0.2S  cm'^  measured  by  ourselves.  The  results 
ploned  in  Fig.3  using  logarithmic  scales  show  two  caracteristic  shoulders  corresponding  to  the 
diflcient  decay  rates. 

1  mm 


cu 

o 

c 

ro 

u5 

I  1/im 

O 


1  ns  1  MS  1  nns 

Time 

Fig.3  Average  distance  xM)  travelled  by  photoexc  .  ed  charge  carriers  as  a  function  of  time  t, 
calculated  fron  Eq.I2  from  our  photocurrent  measurement  at  low  faiences  (^0  pJcm'^)  and 
average  applied  electric  field  of  I  kVcmf 

The  decay  rates  extracted  from  Fig.2  are  consistent  with  the  predictions  that  can  be  made 
unsing  the  parameters  of  Ref.3  for  SUl  and  the  electron-hole  model.  We  take  tj  and  t2  as  the 
respective  hole  and  electron  recombination  times,  associating  the  longest  of  them  with  the 
electrons  which  arc  the  majority  carriers  and  therefore  the  majority  contributors  to  Xafoo)  in  that 
crystal.  Relying  on  the  quantum  efficiency  estimate  of  0.86  from  Ref.3  for  electron  excitation, 

their  mobility  may  be  estimated  to  0.13±0.07  cm^V-ls'*.  The  electron  and  hole  ht  products 
estimated  from  Fig.2  arc  also  in  good  agreement  with  those  determined  holographically  in 
Ref.3.  These  photocurrent  measurements  set  the  proper  experimental  conditions  for  the 
holographic  determination:  The  grating  period  has  to  be  'arge  enough  to  minimize  diffusion  and 

still  within  the  average  travelled  distance  of  100  pm  with  an  applied  field  of  1  kV  .  The  whole 
transport  takes  about  100  psec  which  sets  the  approximate  time  scale  for  the  experiment. 

IV  THE  TIME-OF-FLIGHT  EXPERIMENT  AND  ITS  RESULTS. 

We  experimented  the  holographic  method  using  the  set-up  illustrated  by  figure  4.  The 
light  source  is  the  same  as  in  the  photocurrent  experiment. The  beam  splitter  BSl  provides  two 
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beams  Si  and  S2.  Si  illuminates  a  Ronchi  ruling  R  imaged  on  the  crystal  by  a  pair  of  Fourier 
transforming  lenses  Li  and  L2.  In  the  intermediate  Fourier  plane  F,  an  apenure  selects  two 
lowest  order  spatial  frequency  waves  (  0  and  +1  or  -1)  diffracted  by  R.The  crystal  is  thus 

illuminated  by  a  sinusoidal  distribution  of  intensity  whose  period  is  93  Itm  and  modulation 
depth  m  0  1 .  The  writing  beam  fluence  was  1  pJ/cm^  which  satisfies  the  low  illumination 
condition  (  NA^lO^^m‘3,  ns=2  lO’^cm"^) .  S2  is  used  to  erase  previously  recorded  gratings  by 
uniformly  illuminating  the  crystal  hrom  the  back  .  Two  electronically  controlled  shutters  Shi  and 
Sh2  insure  the  single  pulse  writing  of  a  grating  and  its  erasure.  A  voltage  Vq  is  applied  to  the 
crystal  by  means  of  a  power  supply  which  is  electronically  switched  on  just  before  the  writing 
pulse.The  crystal  is  oriented  for  maximum  diffraction  efficiency  with  the  grating  k  vector  along 
the  (1 10  direction)  and  the  diffraction  efficiency  is  monitored  using  a  10  |lW  He-Ne  laser  beam. 
It  is  linearly  polarized  by-  Pi  so  that  the  polarization  of  the  light,  scattered  by  defects  in  the 
crystal  is  mostly  orthogonal  to  the  diffracted  beam  and  can  be  blocked  by  the  polarizer  P2.  This 
is  made  necessary  by  the  small  grating  period  used  in  this  experiment.  To  minimize  the 
generation  of  charge  carriers  by  the  He-Ne  beam,  an  acousto-optic  modulator  AOM  is  used  to 

turn  the  beam  on  lOOfis  before  the  writing  pulse  from  the  YAG  laser,  and  turn  it  off  during  the 
erasing  phase. 


Fig,4  Experimental  arrangement  for  the  time-of-flight  holographic  mobility  measurement.  ' 

A  carefull  timing  allows  full  decay  of  the  space  charge  field  between  successive  writing  pulse.  A 
P.M.T.  measures  the  time  evolution  of  the  weak  diffracted  signal.  To  minimize  a  space  charge 
field  build-up,  the  voltage  is  connected  about  a  second  before  the  writing  pulse,  while  the  crysud 
is  kept  in  the  dark.  The  noisy  signal  is  averaged  over  about  SO  shots  and  filtered  both 
electronically  at  the  input  of  the  storage  digital  oscilloscope  and  later  numerically.  The  total 
integration  time  of  this  smoothing  process  is  kept  smaller  than  0. 1  times  the  time  at  which  the 
diffraction  peaks  i.e.  the  time  when  Xa(t)  is  half  the  period  of  the  grating. 

Fig.S  shows  the  experimental  time  evolution  of  the  diffraction  efficiency  for  "SUl"  with 

Vq  =  1  kV,  along  with  the  best  fit  of  these  data  to  equation  (11)  with  parameters  ^„„r'andr  ". 
The  real  value  Eq  of  the  electric  field  inside  the  crystal  is  needed  for  the  determination  of  the 
mobility.  It  has  been  frequently  observed  that  Eq  does  not  equal  the  ratio  of  the  applied  voltage 
Vq  to  the  distance  d  between  the  electrodes.  We  measured  it  separately,  using  the  change  in  the 
polarisation  of  the  transmitted  He-Ne  beam  caused  by  the  electrooptic  effect  as  described  in 
Rcf.6.  The  values  measured  for  Eo  were  typically  0.7  Vo/d,  with  the  coefficient  increasing 
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slightly  with  Vq.  To  caracterize  the  uniformity  of  this  field  across  the  crystal,  these 
measurements  were  peifotmed  using  a  beam  diameter  of  0.4  mm  (1/3  of  the  beam  size  used  for 
the  diffraction  experiment).  Eq  was  found  to  be  uniform  within  5  %  both  along  and  normal  to 
the  electric  field  direcdon.  We  also  could  observe  that  the  electric  field  Eq  did  not  change  more 
than  5%  during  the  build-up  of  the  photoiefractive  grating. 


Fig.5  Time  behavior  of  the  diffraction  efficiency  in  SUl  sample  after  a  picosecond  excitation 
pulse.  The  solid  curve  is  a  plot  of  Eq.  11  with  F =8.6  10*  s'^  and  r'=  2.42  lO^s'^  which  gives 
the  best  fit  to  the  data. 


The  mean  electron  drift  velocity  v  =V"lk  .obtained  from  the  fitted  values  of  F'  is  plotted 
on  Fig.6  as  a  funcdon  of  Eq.  The  linear  variadon  confums  the  nodon  of  mobility.  From  this  set 
of  data,  we  obtain  its  value  ^M).24  cm^V'lsec'*.  The  main  uncertainty  remains  here  that  on  the 
electric  field  which  could  be  as  big  as  30%. 
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Fig.6  Mean  electron  drift  velocity  versus  applied  electric  field  as  determined  by  fitting  Eq.l  I  to 
four  curves  as  shown  in  FigJ. 

We  also  nodced  that  contrary  to  Eq.  (10),  the  damping  constant  P  also  varies  linearly 
with  the  electric  field  Eq-  This  can  be  explained  by  a  dispersion  in  the  values  of  the  product  pEo 
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across  the  crystal.  Our  experiments  rule  out  strong  nonuniformities  of  the  electric  field. 
However,  nonuniformities  in  the  mobility  itself  would  easily  result  from  slight  variations  in  the 
environmental  conditions  during  the  growth  of  the  crystal.  This  is  will  be  shown  in  a  future 
publication. 

V  CONCLUSIONS 

In  this  paper  we  demonstrate  the  use  of  an  holographic  method  to  measure  photoexcited  charge 
carriers  mobilities  in  photorefractive  materials.  On  a  well  characterized  n-BSO  crystal,  the 
conductivity  is  seen  to  be  ohmic  over  our  experimental  range  of  electric  field  (200  to  2000  V/cm) 
and  has  a  value  of  0.2410.07  cm2V'*s  '.  We  have  made  two  other  independent  evaluations  of 
the  electron  mobility  in  the  same  sample;  an  electrical  determination  from  the  absolute 
photocurrent,  the  absorption  coefficient  reported  here  and  the  quantum  efficiency  from  Ref.3 
gives  the  value  0.13  ±  0.08  cm^V'U  *.  An  evaluation  based  on  the  measured  diffusion  length 
LD=Ke‘'  =  5.8  ±  2.3  pm  of  Ref.3  and  our  recombination  time  t2=80±  5ps  gives  a  value  0.171 
0.1  cm^V'ls'^  All  three  measurements  give  essentially  the  same  result  which  is  too  small  to  be 
explained  by  collision  limited  free  electron  transport  in  the  conduction  band.  A  model  of 
photoexcited  carriers  hopping  between  thermally  excitable  shallow  traps  as  proposed  in  Ref.7 
would  eventually  explain  those  values.  A  study  of  the  temperature  dependence  of  the  mobility  is 
now  under  progress,  using  the  time-of-flight  technique. 
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We  have  developed  a  time-of-flight  technique  for  measuring  the  mobility  of  photoexcited  charge  car¬ 
riers  in  certain  crystals  exhibiting  the  electro-optic  effect.  We  used  this  holographic  technique  to  find 
that  the  mobility  of  photoexcited  electrons  in  a  previously  well-characterized  sample  of  n-type  Bii2SiO20 
is  0.24±0.07  cm^  V“'  s~',  independent  of  electric  field  in  our  range  of  observation  (»200-2000  V/cm). 

We  also  present  results  of  transient  photocurrent  measurements.  When  used  with  our  absorption  mea¬ 
surement  and  previously  reported  values  of  the  quantum  efliciency  and  mobility-lifetime  product,  they 
give  two  independent  estimates  of  electron  mobility  that  are  consistent  with  our  directly  measured  value. 


I.  INTRODUCTION 

The  mobility  of  photoe.^cited  carriers  in  photorefrac¬ 
tive  insulating  crystals  is  a  key  parameter  in  determining 
the  response  time  of  a  material.  However,  its  reliable 
measurement  has  proven  to  be  more  difficult  than  expect¬ 
ed.  In  nominally  undoped  Bii2SiO20  previous  measure¬ 
ments  using  purely  electrical  experimental  techniques 
have  provided  conflicting  values  ranging  from  5X10“’ 
cm^V“'s“'  to  3  cm^  V“' s“'.'“^  Similarly  conflicting 
values  have  been  reported  for  structurally  similar 
Bii2GeO,2.^’^  The  measurements  of  Refs.  1  and  2  are 
based  on  the  transit  time  of  photoexcited  carriers 
through  a  thin  sample,  but  they  do  not  show  the  usual 
top-hat  profile’  of  photocurrent,  which  is  easy  to  analyze. 
Instead,  shapes  of  current  pulses  lead  to  concepts  like 
“dispersive”  phototransport®  for  which,  strictly  speaking, 
the  mobility  cannot  be  defined.  The  mobility  values  quot¬ 
ed  in  Refs.  3  and  4  are  based  on  the  absolute  magnitude 
of  the  photocurrent  combined  with  measurements  of  the 
optical  absorption  coefficient  and  estimates  of  quantum 
efficiency  (fraction  of  absorbed  photons  creating  photo¬ 
carriers),  thus  giving  only  qualitative  estimates  of  the  mo¬ 
bility.  Purely  electrical  mobility  measurements  also 
suffer  from  complications  caused  by  electrode  potentials, 
surface  currents,  and  nonuniform  electric  fields,  which 
can  even  change  during  the  course  of  the  measurement. 

Optical  measurements  of  the  range  (or  diffusion  length) 
of  photoexcited  charge  carriers^  have  also  been  used  to 
estimate  the  charge-carrier  mobilities  of  Bi,2SiO20  and 
Bii2GeO20-  The  diffusion  length  is  proportional  to  the 
square  root  of  the  mobility  and  the  lifetime  of  the  pho¬ 
toexcited  carriers.  By  measuring  the  lifetime  and  using 
the  diffusion  length  measurement  of  Ref.  7,  Le  Saux,  Lau- 
ney,  and  Brun*  inferred  a  mobility  of  130  cm^  V“’  s~'  for 
carriers  in  Bi|2GeO20.  Jonathan,  Rossignol,  and  Roosen’ 
measured  the  build-up  time  of  the  photorefractive  grating 
excited  by  picosecond  pulses  in  Bi,2SiO20  and  deduced 
the  value  50  cm^  V“ '  s“  ‘  for  the  carrier  mobility.  Anoth¬ 
er  optical  determination  of  the  charge-carrier  mobility 


was  performed  by  Astratov,  Il’inskii,  and  Furman.'® 
They  studied  the  growth  of  electrical  screening  field 
through  the  electro-optic  effect  and  evaluated  values  from 
10“’  to  10“’  cm^V~'s“'  at  temperatures  from  130  to 
200  K,  respectively.  At  higher  temperatures  the  screen¬ 
ing  field  development  was  too  complicated  for  the  deter¬ 
mination  of  mobility." 

Here  we  report  a  direct  holographic  time-of-flight  mea¬ 
surement  of  a  photoexcited  charge-carrier  mobility  in  an 
insulator.'^  We  find  the  electron  mobility  to  be 
0.24±0.07  cm^V“'s“'  in  a  previously  well-charac¬ 
terized  sample  of  Bi,2SiO20-  We  also  report  transient 
photocurrent  measurements  from  which  we  obtain  two 
independent  estimates  of  this  mobility,  (1)  using  the  abso¬ 
lute  value  of  photocurrent,  as  in  Refs.  3  and  4 
(0. 13±0.08  cm^V“'s“'),  and  (2)  using  the  recombina¬ 
tion  time,  as  in  Ref.  8  (0.17±0.1  cm^V“'s“').  These 
two  qualitative  estimates  from  photocurrent  experiments 
are  consistent  with  our  holographic  measurement  of  mo¬ 
bility.  In  the  time-of-flight  mobility  measurement  we  use 
the  spatially  sinusoidal  intensity  pattern  of  overlapping 
picosecond  laser  pulses  to  excite  a  sinusoidal  pattern  of 
electrons  into  the  conduction  band,  where  it  drifts  under 
the  influence  of  a  strong  applied  static  electric  field.  The 
spatially  modulated  space-charge  field  caused  by  the  ex¬ 
cited  electrons  and  the  ions  they  leave  behind  is  initially 
zero,  grows  to  a  maximum  after  electrons  have  moved 
one-half  of  the  period  of  the  intensity  pattern,  and  fall  to 
a  minimum  after  the  full  period.  We  probe  the  space- 
charge  field  with  a  Bragg-matched  cw  laser  through  the 
electro-optic  effect.  The  electron  drift  velocity  and  thus 
the  mobility  are  calculated  from  the  time  required  for 
electrons  to  drift  one  period.  Pauliat  et  al.'^  recently 
used  a  similar  kind  of  method  to  measure  the  mobility  in 
Bi|2GeO20,  however,  with  cw  laser  excitation  and  an  al¬ 
ternating  applied  electric  field. 

We  present  a  simple  theory  for  our  holographic  time- 
of-flight  technique  in  Sec.  II.  Because  estimates  for  mo¬ 
bility  and  lifetime  of  photoexcited  charge  carriers  from 
photocurrent  experiments  were  required  to  set  experi- 
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mental  conditions  suitable  for  the  time-of-flight  mobility 
measurement,  we  describe  the  transient  photocurrent 
measurements  in  Sec.  Ill  before  describing  the  holo¬ 
graphic  mobility  measurement  in  Sec.  IV.  In  Sec.  V  we 
discuss  some  implications  of  the  small  value  of  the  mobil¬ 
ity  (compared  to  drift  mobilities  in  semiconductors), 
which  we  measured. 

II.  THEORY  OF  HOLOGRAPHIC 
TIME-OF-FLIGHT  TECHNIQUE 

We  assume  that  the  photorefractive  crystal  is  il¬ 
luminated  with  two  short  interfering  laser  pulses  that  are 
negligibly  attenuated,  resulting  in  an  energy  flux  distribu¬ 
tion  inside  the  crystal  given  by  (in  units  of  J  cm“^) 

C/(z)=f/o[l  +  Re(me'^)]  .  (1) 

Here  m  is  the  modulation  depth  and  k  is  the  wave  vector 
of  the  interference  pattern.  The  laser  pulse  is  considered 
to  be  so  short  that  the  excited  charge  carriers  may  be  as¬ 
sumed  not  to  move  noticeably  in  space  during  the  pulse. 
We  will  assume  that  after  the  excitation,  the  movement 
of  electrons  is  described  (in  the  dark)  by  the  following 
four  equations.  (Similar  equations  for  holes  may  be  con¬ 
structed  when  needed. ’'*)  First  we  assume  that  the  ion¬ 
ized  donor  density  changes  only  by  the  single 

direct  recombination  process: 

=  (2) 

Here  y  *he  recombination  rate  constant.  Then,  the 
continuity  equation  for  the  density  n  (z,t)  of  photoexcited 
electrons  is 

_  1  3  .  .. . 

of  e  oz 

where  e  is  the  (positive)  charge  of  the  electron.  The 
current  density  ylz.f)  is  assumed  to  consist  only  of  a  drift 
and  a  diffusion  z  component  near  thermal  equilibrium  at 
temperature  T,  so  that  it  can  be  written  in  the  form 


V^(f,z)  =  Vo(f)  +  Re[V,(f)c'*']  , 

M(f,z)  =  noU)  +  Re[«i(f)e'*^]  ,  (6) 

£'(f,z)  =  £o(f)  +  Re[£:,(f)e'^]  . 

Although  we  may  quote  parameters  in  other  units,  the 
formulas  here  will  always  be  assumed  to  be  in  SI  units. 

We  will  assume  that  the  laser-pulse  energy  flux  I/q  is 
small  enough  so  that  Mq  ^^e  proportional  to  Uq. 

Substituting  (6)  in  (2)-(5)  shows  that  this  requires  first 
that  «(,  «N^,  so  that  Np  can  be  replaced  by  ( ==  10** 
cm~^)  in  the  right-hand  side  of  (2).  More  stringently,  it 
requires  the  average  density  of  excited  electrons  to  be 
much  smaller  than  the  saturation  density  defined  by 


so  that  we  only  need  to  keep  terms  in  (3)  that  are  linear  in 
Uq.  In  our  experiments  is  much  less  than 
(  —yN^e/nkgT),  which  was  found  to  be  3X  10*  cm~^  in 
Ref.  14.  This  gives  >  10*^  cm~^  With  these  approxi¬ 
mations  (3)  and  (4)  give 

3n,  kgT 

yN^+k^fj. - --ikfiEQ  ni  ,  (8) 

of  e 

which  can  be  solved  immediately.  With  the  same  as¬ 
sumptions  we  get  an  equation  for  £ , : 


-^  =  7  ^fiEo  +  ikfi—^  j«,  .  (9) 

With  the  solution  for  n,  from  (8)  we  find  a  solution  to  (9): 

£,(f)  =  £„(l-e"''')  ,  (10) 

which  satisfies  the  initial  condition  £](0)==0,  and  in 
which 


r~r'+ir’'=yN^+k  - i^iiEq 


j  =ennE  +  fikgT  -  n  ,  (4) 

oz 

where  /r  is  the  (positive)  electron  mobility  and  kg  is 
Boltzmann's  constant.  Finally,  the  Poisson  equation  for 
the  superposition  £(z,f)  of  CouloTub  fields  arising  from 
the  charged  donor  density  iV/  and  excited  electron  densi¬ 
ty  n  is,  in  SI  (Systeme  International)  units 


where  e  is  the  dielectric  constant  and  the  constant  is 
the  density  of  acceptor  sites  required  to  make  the  crystal 
electrically  neutral  in  the  dark  at  static  equilibrium.  At 
low  excitation  energy  fluxes  the  initial  values  of  Np,  n 
and  £  can  be  written  in  the  same  form  as  the  energy  flux 
of  (1).  Because  we  consider  only  the  case  of  small  modu¬ 
lation  |/n|  «  1,  we  only  need  to  keep  terms  that  are  con¬ 
stant  in  space  or  that  vary  as  e'*'  at  subsequent  times: 


p  Eo+ikkgT/e 

£^  =  -Mjn,(0)-^ - ^ -  .  (12) 

It  is  clear  from  (11)  that  the  imaginary  part  F"  of  the 
damping  constant  does  not  depend  on  any  other  material 
parameters  than  the  mobility  itself.  If  the  experimental 
conditions  are  arranged  to  be  such  that  F"  is  larger  than 
F',  the  temporal  oscillation  in  the  amplitude  £,  of  the 
grating  electric  field  becomes  visible  and  the  mobility  (i 
can  be  determined  from  the  oscillation  period,  i.e.,  the 
time  it  takes  for  an  electron  to  drift  one  grating  period. 
From  (5)  the  oscillation  is  caused  by  the  superposition  of 
the  stationary  ionic  grating  upon  the  grating  of  electrons 
drifting  with  constant  velocity  under  the  influence  of  the 
strong  and  applied  electric  field  Eq,  with  the  two  space- 
charge  fields  canceling  each  other  initially,  and  at  each 
subsequent  time  when  the  gratings  coincide.  The  require¬ 
ment  that  F"  be  clearly  larger  than  F'  corresponds  to  the 
requirement  that  photoexcited  carriers  drift  at  least  one 
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period  before  experiencing  recombination  and  significant 
effects  from  diffusion. 

We  observe  the  time  development  of  the  space-charge 
field  grating  by  diffracting  a  laser  beam  from  it.  The 
diffraction  efficiency  rj  is  proportional  to  the  square  of  the 
amplitude  Ei  of  the  Coulomb  grating,  so  the  time  behav¬ 
ior  of  7)  is,  from  (10), 

V  =  (13) 

where  is  the  value  of  diffraction  efficiency  after  the  os¬ 
cillations  have  damped  out. 

III.  TRANSIENT  PHOTOCURRENT  EXPERIMENTS 

We  studied  photoexcited  carriers  in  a  5X5X5  mm^  n- 
type  Bii2SiO20  single  crystal  grown  by  Sumitomo  and 
designated  as  crystal  SUl  in  Ref.  14,  where  it  was  charac¬ 
terized  by  fitting  quasi-cw  two-beam  coupling  and  grating 
erasure  data  to  the  electron-hole  competition  model. 
These  experiments  did  not  determine  the  mobility  fi  and 
the  recombination  time  r  separately,  but  only  their  prod¬ 
uct.  Therefore  we  first  made  the  following  photocurrent 
measurements  to  obtain  estimates  of  fi  and  r  separately, 
and  to  predict  a  useful  grating  spacing  for  the  time-of- 
flight  measurement. 

We  studied  the  decay  of  the  photocurrent  induced  in 
the  crystal  by  a  uniformly  illuminating  30-ps  laser  pulse 
when  an  electric  field  ( =  1  kV/cm)  was  applied  to  the 
crystal,  which  is  coated  with  two  conducting  electrodes 
separated  by  d.  At  low  fluences  (  5  30  fiJ/cm^)  the  decay 
of  the  photocurrent  i  ( t)  can  be  fit  by  two  exponential  de¬ 
cays  over  several  orders  of  magnitude  in  time  and  level  of 
signal  (as  reported  in  Ref.  3).  The  time  constants  for  the 
exponential  decays  are  r|=26±3  ns  and  /2~80±5^s. 
The  faster  component  displaced  about  1%  of  the  charge 
of  the  slower  one.  A  convenient  way  to  present  the  pho¬ 
tocurrent  data  is  to  express  it  in  the  form  of  the  average 
magnitude  Xglt)  of  the  distance  traveled  by  a  photoexcit¬ 
ed  charge  carrier: 

x,it)  =  N-'  i\x„it)\^^fW)df  (14) 

and  plot  the  curves  using  logarithmic  scales.  Here  x„(t) 
is  the  distance  traveled  by  the  nth  carrier  in  the  direction 
of  the  electric  field,  and  N  is  the  total  number  of  mobile 
carriers  generated  by  the  laser  pulse.  We  estimate  N  by 
applying  the  observation  of  Ref.  14  that  nearly  every  ab¬ 
sorbed  photon  excites  a  carrier  to  our  measured  absorp¬ 
tion  coefficient  of  0.74±0.25  cm~'.  Taking  into  account 
the  20%  Fresnel  reflection  at  each  surface  of  the  crystal, 
we  find  the  Intx^,)  vs  ln(t)  curve  shown  in  Fig.  1  with  the 
two  characteristic  shoulders  corresponding  to  the  two 
different  exponential  decays.  I  he  parameters  of  the  two 
exponentials  are  consistent  with  the  predictions  that  one 
would  make  assuming  the  electron-hole  model  and  the 
parameters  of  Ref.  14  for  SUl,  if  we  take,  in  addition,  the 
hole  and  electron  recombination  times  to  be,  respectively, 
the  t|  and  <2  of  F'g-  1-  However,  any  reasonable  model 
would  have  to  associate  the  longer  t2  with  the  electrons 
that  are  clearly  the  majority  carriers  for  the  photorefrac- 
tive  effect,  and  the  majority  contributors  for  x„(oc)  in 
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FIG.  1.  Average  distance  x„(/)  traveled  by  photoexcited 
charge  carriers  as  a  function  of  time  /,  calculated  by  Eq.  (14) 
from  our  phctocurrent  measurements  at  low  fluences 
( 1 30  fiJ/cm')  and  at  average  applied  electric  field  of  1  kV/cm. 

Fig.  1.  We  see  that  the  majority  carriers  (electrons)  travel 
about  100  fj.m,  before  recombination  when  the  applied 
electric  field  is  of  the  order  of  1  kV/cm.  Therefore  we 
adjust  the  grating  period  to  be  about  100  fim  in  the  ex¬ 
periments  of  the  next  section. 

We  can  determine  a  preliminary  value  for  the  electron 
mobility  [fj.^  =Xa(t)AEQt)]  from  Fig.  1.  Using  values  for 
Xa  and  t  just  before  the  second  shoulder  and  taking  into 
account  the  quantum  efficiency  of  0.86  for  electron  exci¬ 
tation'^  we  get  fi=0. 13±0.08  cm^  V'  s“'.  As  was  men¬ 
tioned  in  the  Introduction,  another,  independent  estimate 
of  mobility  uses  the  recombination  time  /2  =  80±5^s 
measured  here  and  diffusion  length  A'J"'  =  5.8±2.3  ^m 
reported  in  Ref.  14  for  electrons  in  the  same  sample. 
This  gives  a  value  /i  =  0. 17±0. 1  cm"'  V~'  s~'. 

IV.  HOLOGRAPHIC  TIME-OF-FLIGHT 
MOBILITY  EXPERIMENT 

Our  holographic  time-of-flight  mobility  measurement 
employs  the  experimental  arrangement  illustrated  in  Fig. 
2.  The  interference  pattern  is  written  by  30  ps 
frequency -doubled  pulses  (at  5  Hz)  from  a  Nd/YAG  (yt¬ 
trium  aluminum  garnet)  laser  at  the  wavelength  of  532 
nm.  The  beam  expander  (BE)  creates  uniform  illumina¬ 
tion  of  the  sample.  The  beam  splitter  BSl  provides  two 
beams  5 1  and  Si-  Beam  S,  illuminates  a  Ronchi  ruling  R 
which  is  imaged  on  the  Bi,2SiO20  sample  by  a  pair  of 
Fourier-transforming  lenses  L,  and  Lj-  ^he  intermedi¬ 
ate  Fourier  plane  an  aperture  A  selects  the  two  lowest- 
order  components  (0  and  —  1  or  1)  diffracted  by  R.  The 
crystal  is  thus  illuminated  by  a  sinusoidal  distribution  of 
intensity  whose  period  is  93  ^m  and  modulation  depth 
|m|  is  near  1.  The  writing  pulse  energy  flux  was  1 
ftJ/cm^  (nQ=1.6X10'^  cm“^).  Beam  Sj  is  used  to  erase 
the  previously  written  gratings  by  uniformly  illuminating 
the  crystal  from  the  back.  Two  electronically  controlled 
shutters  Sh,  and  Shj  ensure  single-pulse  writing  of  the 
grating  and  its  erasure.  A  voltage  Kg  is  applied  to  the 
crystal  by  a  power  supply  which  is  electronically 
switched  on  just  before  the  writing  pulse  to  minimize 
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FIG.  2.  Experimental  arrangement  for  time-of-flight  mea¬ 
surement  of  mobility  of  electrons  in  sample  Bi,2SiO20. 

space-charge  field  buildup.  The  crystal  is  oriented  so  that 
the  grating  wave  vector  is  along  its  (110)  crystallographic 
axis.  The  diffraction  efficiency  is  monitored  using  a  10- 
/iW  He-Ne  laser  beam  that  is  linearly  polarized  by  so 
that  the  polarization  of  the  light  scattered  by  defects  in 
the  crystal  is  largely  orthogonal  to  that  of  the  diffracted 
beam  and  can  be  blocked  by  the  polarizer  ^2-  ^.n 
acousto-optic  modulator  (AOM)  is  used  to  turn  the  He- 
Ne  beam  on  synchronously  with  the  writing  pulse  from 
the  Nd/YAG  laser,  and  turn  it  off  during  the  erasure 
phase.  A  careful  timing  allows  full  decay  of  the  space- 
charge  field  between  successive  writing  pulses  and  mini¬ 
mizes  the  effect  of  the  photocurrents  produced  by  the 
He-Ne  beam,  which  is  switched  on  only  200  us  before  the 
writing  pulses.  A  photomultiplier  tube  (PMT)  measures 
the  time  evolution  of  the  diffracted  signal.  The  noisy  sig¬ 
nal  is  averaged  over  about  50  shots  and  filtered  both  elec¬ 
tronically  at  the  input  of  the  storage  digital  oscilloscope 
and  later  numerically.  The  total  integration  time  of  this 
smoothing  process  is  kept  smaller  than  0. 1  times  the  time 
interval  between  the  laser  pulse  and  the  diffraction  max¬ 
imum. 

Figure  3  shows  the  experimental  time  evolution  of  the 
diffraction  efficiency  for  SUl  with  the  voltage  Fo=l  kV 
applied  across  the  crystal,  along  with  the  best  fit  of  this 
data  to  (13)  having  T',  and  F"  as  fitting  parameters. 
To  be  able  to  evaluate  the  mobility  from  (13)  we  must 
know  the  value  of  the  applied  electric  field  £q  inside  the 
crystal.  It  has  been  frequently  observed  that  Eq  inside 
photorefractive  materials  does  not  generally  equal  f^o/d, 
where  d  is  the  distance  between  the  electrodes.”  We 
measured  the  electric  field  Eq  in  a  separate  experiment 


0  20  40  60 


t  (/iS) 

FIG.  3.  Time  behavior  of  diffraction  efficiency  in  Bii2SiO20 
sample  after  picosecond  photoexcitation  pulse.  The  solid  curve 
is  a  plot  of  Eq.  (13)  with  r'  =  8.6X  10^  s“‘  and  r"  =  2.42X  10’ 
s~',  which  gives  the  least  mean-squares  error  when  compared  to 
data. 

using  the  change  in  the  polarization  of  the  transmitted 
He-Ne  beam  caused  by  the  electro-optic  effect  and  using 
the  analysis  of  Ref.  15.  The  values  measured  for  Eq  were 
typically  0.7  Vq  /d,  with  the  coefficient  increasing  slightly 
with  Vq.  To  characterize  the  uniformity  of  Eq  across  the 
crystal  these  measurements  were  done  with  a  beam  diam¬ 
eter  of  0.4  mm,  one-third  of  the  beam  size  used  for  the 
diffraction  experiment.  The  field  Eq  was  found  to  be  uni¬ 
form  within  5%  both  along  and  normal  to  the  electric- 
field  direction.  We  also  could  observe  that  Eq  did  not 
change  more  than  5%  during  the  buildup  of  the  pho¬ 
torefractive  grating. 

The  imaginary  part  T"  of  the  parameter  in  (11)  divided 
by  k  is  the  mean  electron  drift  velocity  v.  A  plot  of  the 
drift  velocity  v,  obtained  from  best-fit  T"  values  as  a 
function  of  Eq  is  shown  in  Fig.  4,  demonstrating  linear 
variation  and  confirming  that  the  photoconduction  is 
Ohmic.  From  this  set  of  data  we  obtain  a  mobility 
fi=0.2A  cm^V“'s“'  with  the  main  uncertainty  arising 
from  the  uncertainty  in  the  average  electric  field,  which 
we  estimate  could  be  as  far  off  from  our  electro-optic 
measurement  as  is  Kg/d  (  ==±30%).  We  also  found  that, 
contrary  to  (11),  the  best-fitting  T'  also  varies  linearly 
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FIG.  4.  Mean  electron  drift  velocity  v  vs  applied  electrc  field 
Eq  as  determined  by  fitting  Eq.  (13)  to  four  curves,  such  as 
shown  in  Fig.  3.  iJ  =  r" /k. 
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with  the  electric  field  Eq  for  points  of  Fig.  4.  We  have 
shown  that  this  is  consistent  with  a  random  variation  of 
mobility  ( ±40%)  across  the  crystal  volume.’^ 

V.  DISCUSSION 


The  electron  mean  free  path,  calculated  from  the  stan¬ 
dard  expression  [^(2Agrm)’''^/e]  and  our  measured 
/I  =  0.24  cm^  V'  '  s“',  is  less  than  the  typical  separation 
of  atoms  in  the  crystal  lattice.  This  indicates  that  the 
model  of  collision-limited  free-electron  transport  in  the 
conduction  band  does  not  apply  to  the  photoexcited  elec¬ 
trons  that  we  are  observing.  However,  our  measurements 
demonstrate  that  the  average  drift  velocity  depends 
linearly  on  the  applied  electric  field  (Fig.  4)  and  therefore 
justifies  the  concept  of  mobility.  Our  low  mobility  value 
can  be  explained  as  a  trap-limited  mobility  Pho¬ 

toexcited  electrons  do  not  move  in  the  conduction  band 
continuously.  They  fall  occasionally  to  an  energy  level 
from  which  they  can  be  thermally  excited.  The  trap- 
limited  mobility  is  related  to  the  actual  conduction-band 
mobility  by  the  expression 


where  t,  is  the  average  time  spent  by  a  conduction  elec¬ 
tron  between  trapping  events  and  r,  is  the  average  time 
for  a  trapped  electron  before  its  release.  Shallow-trap 
models  have  already  been  used  in  analyses  of  charge 
transport  in  piiot'^refractive  materials.^”'’  '*’  ’’ 

The  wide  range  of  mobility  values  reported  in  the 
literature  may  also  be  an  indication  of  trap-limited  trans¬ 


port  processes;  small  differences  (greater  than  kg  7)  in 
the  depths  of  traps  can  indeed  lead  to  an  order-of- 
magnitude  difference  in  thermal  excitation  rates  and  thus 
in  trap-limited  mobilities.  By  their  nature,  shallow-trap 
events  are  always  dominated  by  a  very  narrow  energy  in¬ 
terval  in  the  gap.'* 

In  conclusion,  we  demonstrate  here  how  a  holographic 
method  can  be  used  to  measure  mobilities  of  photoexcit¬ 
ed  charge  carriers  in  photorefractive  materials.  Using 
this  method  on  a  well-characterized  n-type  Bi,2SiO20 
crystal  we  found  a  mobility  of  0  24±0.07  cm^  V” '  s“ '  for 
the  electrons.  The  conductivity  caused  by  these  photoex¬ 
cited  electrons  was  found  to  be  Ohmic  over  our  experi¬ 
mental  range  of  200-2000  V/cm  in  the  electric  field.  We 
also  estimated  electron  mobility  values  indirectly  by  two 
independent  ways  from  transient  photocurrent  measure¬ 
ments  reported  here  and  from  previous  holographic  char¬ 
acterization  of  the  same  sample.  We  found  good  agree¬ 
ment  among  all  three  mobility  values. 
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